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Abstract 

Using the large acceptance apparatus FOPI, we study central collisions in the re- 
actions (energies in A GeV are given in parentheses): ^'^Ca+'^'^Ca (0.4, 0.6, 0.8, 
1.0, 1.5, 1.93), ^^Ni+s^Ni (0.15, 0.25, 0.4), ^'^Ru+^'^Ru (0.4, 1.0, 1.5), 962r+962r 
(0.4, 1.0, 1.5), I29xe+Csl (0.15, 0.25, 0.4), ^^'^Au+^^'^Au (0.09, 0.12, 0.15, 0.25, 0.4, 
0.6, 0.8, 1.0, 1.2, 1.5). The observables include cluster multiplicities, longitudinal 
and transverse rapidity distributions and stopping, and radial flow. The data are 
compared to earlier data where possible and to transport model simulations. 

Key words: heavy ions, rapidity, stopping, viscosity, radial flow, cluster 
production, nucleosynthesis, nuclear equation of state, isospin 
PACS: 25.75.-q, 25.75.Dw, 25.75.Ld 



1 Introduction 



In the past twenty five years studies of nucleus-nucleus collisions over a very broad range 
of energies have been performed. The goal was, and still is, to learn more about the prop- 
erties of hot and dense nuclear matter, namely its equation of state (EOS) and transport 
coefficients such as viscosity and thermal conductivity. In the energy range (0.1 to 2 A GeV) 
of the Schwerionen Synchrotron (SIS) accelerator located at GSI, Darmstadt, one expects 
from simple-minded hydrodynamic considerations that densities twice to three times sat- 
uration density are reached. Another important characteristic of this energy range is that 
the compressed overlap zone created in the collision expands at a rate comparable to the 
rate at which the accelerated nuclei interpenetrate, a phenomenon that provides a very 
useful clock: the passage time 2R/{pcm/'mN) where R is the nuclear radius, pcm and niN 
are the incoming nucleon's cm. momentum and mass. In a historic paper [1] a strong 
ejection of matter in the transverse direction to the beam was predicted. In the present 
work on very central collisions, this transverse direction, in particular in comparison to 
the longitudinal (beam) direction, will play a crucial role. 



The discovery p|3] of collective flow seemed to confirm that fiuid dynamics was the proper 
language to theoretically accompany the experimental efforts. However, soon after it be- 
came clear that a quantitative extraction of fundamental properties of nuclear matter 
required the development of microscopic transport theories [1] that did not have to rely 
on the local equilibrium postulate needed to apply the (one-fluid) hydrodynamic approach 
[5115] . Immediately it also became then clear that the static EOS was not the only unknown 
as momentum dependent potentials [7118] and in-medium modiflcations of hadron masses 
[9j and cross sections [lOj came into play. Despite of these increased difficulties the past 
decade has witnessed reasonably successful attempts to come to a first-generation of nu- 
clear EOS using experimental observations in heavy ion collisions as constraints 
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One approach, [TT], choose to base its conclusions on azimuthally asymmetric flow mea- 
surements for the Au+Au system (reviewed in p3fT5] ) done at the BEVALAC and AGS 
accelerators (USA), while the other approach [T21IT6] was based on a comparison of kaon 
production in C+C and Au+Au collisions using primarily data obtained at SIS by the 
KaoS Collaboration ([IT], for a review see [T81IT9] ). 

Despite this encouraging progress a critical look at the present situation reveals that more 
work is necessary. Fundamental conclusions based on one particular observable, such as 
the K^ production varying the system size [IT] or centrality [T6], are not sufficiently con- 
vincing to settle the problem. At SIS energies kaons are a rare probe: besides the condition 
that their production in a dense medium must be under sufficient theoretical control, one 
must also assume that the change of environment, i.e. the bulk matter produced, is cor- 
rectly described. This involves a complete mastery of the degree of stopping reached in the 
collisions, since the idea |20] behind the kaon observable is that it is (in the subthreshold 
regime) sensitive to the achieved density. Further, the possible influence of isospin depen- 
dences has to be accounted for: the switch from the iso-symmetric system ^^C+^^C to the 
asymmetric ^^^Au+^^^Au is expected to produce a shift of the ratio of the isospin pair 
K^ / K^ , but the K^ production was not measured. 

On the other hand, EOS determinations based on reproducing flow data are also not 
straightforward as has become clear from a number of earlier works by our Collaboration 
[2T]I22]I28]IM1I25]I26]I271I28]^ and the EoS Collaboration [30]. Estimates using Rankine- 
Hugoniot-Taub shock equations [131I31J . (that ignore geometrical constraints and non- 
equilibrium, as well as momentum-dependent phenomena) , show that the pressures gener- 
ated in the collisions, which are at the origin of the measured flow observables, differ only 
relatively modestly when the assumed stiffness of the 'cold' (zero temperature) EOS is 
varied Lil, leading already in 1981 the authors of ref. [31] to conclude that high accuracy of 
flow and entropy measurements would he needed to deduce information on the cold EOS. In 
this connection it is interesting to note that a promising new observable has been proposed 
[28] suggesting to use azimuthal dependences of the kinetic energies, rather than the yields 
of emitted particles. A disadvantage of the azimuthally asymmetric flow signal, is that it 
converges to zero for the most central collisions where compression is maximal. 

It is also a somewhat unsatisfactory situation, that in the past no clean isolation of the 
isospin dependence of the EOS has been made: most of the conclusions based on heavy 
ion data involve the Au+Au system, which is neither isospin symmetric nor close to the 
composition of neutron matter (or neutron stars). 

In the present work, devoted to hulk observations in very central collisions, we set up a 
large data base of additional challenging constraints for microscopic transport simulations 
that claim to extract basic nuclear matter properties from heavy ion reactions. We shall 
show that sensitivity to the EOS is not limited to rare probes or non-central collisions. 

One observable, the ratio of variances of transverse and longitudinal rapidities [32], that 
we shall loosely speaking call 'stopping', will catch our particular attention. Developing 
further on the original idea [T], we shall show evidence that this ratio is influenced by 



it is primarily the density that is modified 



the ratio of the speed of expansion and the speed at which the incoming material is 
progressing and in this sense the 'clock' is not limited to non-central collisions. A first 
attempt to use simultaneously the 'stopping' observable and the directed flow to constrain 
the EOS was published in [33]. Another signal of the EOS that we see is the degree of 
cooling following the decompression of the system: correlations between stopping, radial 
flow and clusterization serve as evidence that at least a partial memory of the original 
compression exists [M]. Cooling (reabsorption for produced particles, clusterization for 
nucleonic matter) can be more efficient if it was preceded by more compression (softer 
EOS). This has different, subtle, consequences for pion, kaon and nucleonic cluster yields 
at freeze out. 

It is well known by now that this kind of physics has overlap with some areas of astro- 
physics. The stiffness of the EOS limits the density in the center of neutron stars as it 
does limit the maximum densities achievable in heavy ion reactions. The synthesis of nuclei 
during the expansion-cooling processes in the Universe and more locally of supernovae is 
partially imitated in the expansion-cooling of the heavy-ion systems we observe. 

With this exciting context in mind we have studied twenty five system-energies varying 
the incident energy by a factor twenty and choosing three, and for some incident energies 
more, different system sizes and compositions in an effort to elucidate the finite size and 
isospin effects on the observables: it is well known that nuclear ground state masses are 
heavily affected by direct (surface) and indirect (Coulomb and neutron-proton asymmetry) 
contributions due to finite size: the 'correction' to the bulk energy is by a factor two. Nuclei 
are therefore truly 'small' objects. In the case of nuclear reactions with incident energies 
not larger than the rest masses an added finite dimension comes to the nuclear 'surface': 
the finite time duration. There is no reason to expect then, that finite size effects are 
smaller in reaction physics than for ground state masses. Data analyses, in particular 
thermal model analyses, that ignore this fact are likely to be misleading. 

Relative to other earlier experimental studies we shall make a substantial effort to elimi- 
nate, or at least clarify, apparatus specific effects on the observables. First, the choice of the 
centrality selection is done in such a way that a precise matching of centralities between 
experiment and simulation is possible: this is important, as some bulk observables (yields, 
stopping) turn out to be very sensitive to centrality for good physical reasons. Second, we 
try to make superfluous the use of apparatus fllters when comparing data from experiment 
and simulation. A large effort is devoted to reconstruct the full An data wherever this is 
possible with a reasonably low uncertainty. The large acceptance of our setup, the FOPI 
apparatus [35|36] . makes this possible. 

We start our paper with a more technical part, a brief description of the apparatus, section 
121 followed by some frequently used deflnitions and the centrality selection, section [31 and a 
detailed account of our Air reconstruction method, section |H The description of results will 
then start with a presentation in section [5] of a systematics of longitudinal and transverse 
rapidity distributions for various identifled ejectiles. We proceed with deduced information 
on the three stages of the reaction: compression, or stopping, section El expansion, or 
radial flow, section [7] and freeze-out, or chemistry, section |H1 We will use the transport 
code IQMD [37] both for technical checks and for an assessment of the basic information 
that can be extracted from our data. References to some of the copious earlier literature 



will be presented wherever relevant in the various sections. We provide a summary (besides 
the abstract and the definition section) for the quick reader. 

The spirit of the present work is similar to the one on pion emission at SIS energies [38] and 
another paper in preparation on azimuthally asymmetric fiow: provide a large systematic 
data base (for the SIS energy range) that hopefully can be used for a reassessment of a 
second generation of nuclear EOS and nuclear transport properties, accompanied by the 
hope that the various existing or still developed transport codes will eventually converge 
to common conclusions. 



2 Apparatus and data analysis 

The experiments were performed at the heavy ion accelerator SIS of GSI/Darmstadt using 
the large acceptance FOPI detector [35l|36j . A total of 25 system-energies are analysed for 
this work (energies per nucleon, E/ A, in GeV are given in parentheses): '^"Ca+^'^Ca (0.4, 
0.,6 0.8, 1.0, 1.5, 1.93), s^Ni+s^Ni (0.15, 0.25), ^^Sxe+Csl (0.15, 0.25), ^^Ru+^^Ru (0.4, 
1.0, 1.5), 962r+96zr (0.4, 1.5), ^^^Au+^^^Au (0.09, 0.12, 0.15, 0.25, 0.4, 0.6, 0.8, 1.0, 1.2, 
1.5). 

Two somewhat different setups were used for the low energy data {E/A < 0.4 GeV) and the 
high energy data. For the latter, particle tracking and energy loss determination were done 
using two drift chambers, the CDC (covering polar angles between 35° and 135°) and the 
Helitron (9° — 26°), both located inside a superconducting solenoid operated at a magnetic 
field of 0.6 T. A set of scintillator arrays. Plastic Wall (7° — 30°), Zero Degree Detector 
(1.2° -7°), and Barrel (42° -120°), allowed us to measure the time of fiight and, below 30°, 
also the energy loss. The velocity resolution below 30° was (0.5 — 1.5)%, the momentum 
resolution in the CDC was (4 — 12)% for momenta of 0.5 to 2 GeV/c, respectively. Use of 
CDC and Helitron allowed the identification of pions, as well as good isotope separation 
for hydrogen and helium clusters in a large part of momentum space. Heavier clusters are 
separated by nuclear charge. More features of the experimental method, some of them 
specific to the CDC, have been described in ref. [39] . 

In the low energy run the Helitron drift chamber was not yet available and the Barrel 
surrounding the CDC covered only 1/8 of the full 27r azimuth. On the other hand a set 
of gas ionization chambers (Parabola) installed between the solenoid magnet (enclosing 
the CDC) and the Plastic Wall was operated that extended to lower velocities and higher 
charges the range of charge identified highly ionizing fragments and covered approximately 
the same angular range as the Plastic Wall. The Parabola together with the Plastic Wall 
represent what we used to call the early PHASE I of FOPI the technical details of which 
can be found in ref. [35]. Also in the low energy run, the CDC was operated in a 'split mode' 
in an effort to extend the dynamical range of measured energy losses along the particle 
tracks: of the 60 potential wires the inner 30 wires were held under a lower voltage, see ref. 
[23] for more details. We note that the combination of the two setups allowed us to cover 
a rather large range of incident E/A: there is approximately a factor of twenty between 
the lowest and the highest energy. 



3 Some definitions, centrality selection 



Choosing the cm. as reference frame, orienting the z-axis in the beam direction, and 
ignoring deviations from axial symmetry the two remaining dimensions are characterized 
by the longitudinal rapidity y = 1/^, given by exp{2y) = {l + f3z)/{l — f3z) and the transverse 
(spatial) component t of the four-velocity u, given by Ut = f3t'y. The 3-vector (3 is the 
velocity in units of the light velocity and 7 = 1/a/1 — f3'^. In order to be able to compare 
longitudinal and transversal degrees of freedom on a common basis, we shall also use the 
transverse rapidity, y^, which is defined by replacing P^ by (3^ in the expression for the 
longitudinal rapidity. The x-axis is laboratory fixed and hence randomly oriented relative 
to the reaction plane, i.e. we average over deviations from axial symmetry. The transverse 
rapidities y^ (or yy) should not be confused with yt which is defined by replacing (3z by 

For thermally equilibrated systems /3j = y/2f3x and the local rapidity distributions dN/dy^ 
and dN/dyy (rather than dN/dyt) should have the same shape and height than the usual 
longitudinal rapidity distribution dN/dyz, where we will omit the subscript z when no 
confusion is likely. Throughout we use scaled units yo = y/yp and Uto = Ut/up, with 
''J'P = Pplp) the index p referring to the incident projectile in the cm.. In these units the 
initial target-projectile rapidity gap always extends from j/q = ~1 to y^ = 1. Besides 
the 'scaled transverse rapidity', y^o, distributions dN/dy^o we will also present data on 
'constrained scaled transverse rapidity' distributions dN/dyxmo which are obtained using 
a midrapidity-cut \yzo\ < 0.5 on the longitudinal rapidities yzo- In this work the shapes 
of the various rapidity distributions will often be characterized by their variances which 
we term varz, varx, and varxm for the longitudinal, transverse and constrained transverse 
distributions, respectively, adding a suffix zero if scaled units are used. In section 15.31 we 
will introduce the ratio varxz = varx /varz which we propose as a measure of stopping to 
be discussed in detail in section El This ratio, naturally, is scale invariant. 

Collision centrality selection was obtained by binning distributions of either the detected 
charge particle multiplicity, MUL, or the ratio of total transverse to longitudinal kinetic 
energies in the center-of-mass (cm.) system, ERAT. To the degree where ERAT is evi- 
dently less influenced by clusterization than MUL, there is some advantage of using ERAT 
when trying to simulate the selections with transport codes that do not reproduce well the 
measured multiplicities. As will be clear from our data, a careful matching of centralities 
between data and simulations is important for correct quantitative conclusions. ERAT is 
related to the stopping observable varxz defined earlier. To avoid autocorrelations we have 
always excluded the particle of interest, for which we build up spectra, from the definition 
of ERAT. 

We estimate the impact parameter h from the measured differential cross sections for the 
ERAT or the multiplicity distributions, using a geometrical sharp-cut approximation. The 
ERAT selections show better impact parameter resolution for the most central collisions 
than the multiplicity selections. In the present work we will limit ourselves to ERAT 
selected data. More detailed discussions of the centrality selection methods used here can 
be found in refs. [33|l0] . In the sequel, rather than using cryptic names for the various 



centrality intervals, we shall characterize the centrality by the interval of the scaled impact 
parameter bo defined by Bq = b/b^ax- We take b^ax = 1.15(y4p + Ar/ ) fm giving values 
known |41] to describe well effective sharp radii of nuclei with mass higher than 40u. This 
scaling is useful when comparing systems of different size. In this work we limit ourselves 
to four centrality bins: 60 < 0.15, 60 < 0.25, 0.25 < 60 < 0.45 and 0.45 < feo < 0.55. Most 
of the data presented here concern the most central bin. In this most central bin typical 
accepted event numbers vary between (45 — 90) x 10^ for all but the Ca+Ca data where we 
registered (20 — 35) x 10'^ most central events. The total number of registered events was 
approximately ten times higher and was limited by a multiplicity filter. Some minimum 
bias events were registered at a lower rate. 

When comparing to IQMD simulations we also take the ERAT observable for binning and 
carefully match the so defined centralities to those of the experiment. Apparatus effects 
tend to modify the observed shape of the ERAT distributions. However, the influence of 
a realistic apparatus fllter on the ERAT binning quality was found to be very small and 
was therefore not critical as long as the corresponding cross sections remained matched. 

The IQMD version we use in the present work largely corresponds to the description given 
in Ref. [37]. More speciflcally, the following 'standard' parameters were used throughout: 
L = 8.66/m^ (wave packet width parameter), t = 200 fm (total integration time), K = 
200, 380 MeV (compressibility of the momentum dependent soft, resp. stiff EoS, Egym = 
25p/po MeV (symmetry energy, with po the saturation value of the nuclear density p). 
The versions with K = 200, resp. 380 MeV are called IQMD-SM, resp. IQMD-HM. The 
clusterization was determined from a separate routine using the minimum spanning tree 
method in conflguration space with a clusterization radius Re = 3 fm. 



4 An reconstruction 



It is of high interest to know the event topology in full Att coverage. Thus, Air yields are 
needed for convincing estimates of chemical freeze-out characteristics, chemical temper- 
atures, entropy generation etc. Also stopping characteristics to be extensively discussed 
later, require a large coverage of phase space. Since measured momentum space distribu- 
tions do not cover the complete An phase space they must be complemented by interpo- 
lations and extrapolations. As this is a rather important feature of our data analysis, we 
devote a whole section to this undertaking. 

Since this study is limited to symmetric systems, we require reflection symmetry in the 
center of momentum {cm.). 

To start the data treatment, we fllter the data to eliminate regions of distorted measure- 
ments, such as edge effects. One advantage of flltering experimental data, rather than 
just the simulated events, is that well deflned sharp limits are introduced which can then 
be applied in exactly the same form to simulations (see later). Examples of these sharp 
flltered data are given in Figs. 1,2,4,5 and 6. A second advantage of these sharp fllters, 
that are more restrictive than just the nominal geometries of the various sub-detectors, is 
that within the remaining acceptances detector efficiencies are close to 100%. For exam- 



pie, studies of measured correlated adjacent detector hit probabilities in the Plastic Wall 
showed a worst case (Au+Au at 1.5 A GeV) double hit probability of 8%, lowering the ap- 
parent total multiplicity, which was nearly compensated by a cross-talk between detector 
units of 6%, raising the apparent multiplicity. This cross talk was mainly due to imperfect 
alignment of the beam axis with the apparatus. The rationale for dropping further going 
detailed efficiency considerations within the filtered acceptances was that the charge bal- 
ances, after in reconstruction, were found to be accurate to typically 5% (see the tables 
in the appendix) for all the 25 system-energies studied, i.e. for a statistically significant 
sampling encompassing a large variation in energy and multiplicity. In our analysis the 
Plastic Wall and the CDC are treated as 'master' detectors required for the registration of 
an ejectile. The Helitron efficiency was taken care of by matching its signals for H isotopes 
to Z=l hits in the Plastic Wall after subtracting from the Plastic Wall hits the estimated 
pion contributions known from our earlier study [38] • At incident energies beyond lA GeV, 
where the Plastic Wall showed a large background for Z > 2, a matching of Z with the 
Helitron was required (in addition to track matching) using the energy loss signals and 
assuming for the Helitron the same efficiency as for Z=l. Other assumptions were found to 
be in conflict with the global charge balances. The Barrel, as a slave detector to the CDC, 
was used to improve isotope identification (for H and He) whenever its signal matched the 
CDC tracking. The mixing of adjacent particles is estimated to be less or equal to 10%, 
except for tritons and ^He in the low energy run where it could be up to 20%. 

In the following step each measured phase-space cell dyo x duto and its local surrounding 
Nyu — 1 cells, with Ny^ = {2ny + l)(2n„ + 1); is least squares fitted using the ansatz 

1 d'N 

exp[/(yo,Mto)J 



Mjo dutodyo 

where f{yo, Wto) = «2Z/o + ^s^^g + d\yo\uto + ai\yo\ + biUto + cq is a five parameter function. 

This procedure smoothens out statistical errors and allows subsequently a well defined 
extension to gaps in the data. Within errors the smoothened representation of the data 
follows the topology of the original data: deviations of 5% are caused by local distortions 
of the apparatus response, revealing typical systematic uncertainties. 

The technical parameters of the procedure were chosen to be dyQ = 0.1, duto = 0.04, Uy = 2, 
riu = 5, i.e. the local fitting domain consisted of a maximum of 55 cells. These choices 
were governed by the available statistics and the need to follow the measured topology 
within statistical and systematic errors. Variations of these parameters were investigated 
to determine the systematic errors of the procedure. 

The smoothened data are well suited for the final step: the extrapolation to zero transverse 
momenta. The low pt extrapolation was extensively checked by event simulations using 
both thermal models or transport codes. For fragments separated only by charge the phase 
space coverage of FOPI is rather good and the extension to An is unproblematic. When 
reconstructing the isotope separated distributions we have used the constraint that the 
sum of isotopes of Z=l and 2 should be consistent in every phase space cell dyo x duto 
with the total value for Z=l, resp. Z=2 fragments (henceforth dubbed 'sum-of-isotopes 
check'). This ensures at low pt that Coulomb effects are approximately respected. 



Finally, a word on systematic errors: Systematic errors are dominated by forward-backward 
inconsistencies and extrapolation uncertainties. The estimates for the latter were guided by 
the accuracy of the total charge balances and the isotope balances. The global systematic 
errors on particle yields are given in the appendix. The small straggling of differential 
data due to generally good statistics and the smoothening fits suggest point to point 
errors smaller than the global errors. 

In the following we show a few illustrative examples of the Aw reconstructions. 
4-1 Thermal model reconstruction 



The above ansatz in terms of exponential functions to limited two-dimensional phase space 
cells is inspired by the thermal model, but does allow for deviations from it. Of course 
in case the thermal model would represent the data well, we would recover the thermal 
model parameterization (i.e. the 'temperature'). This is shown in our first example below 
(Fig. H]). The sharp filter used in this case is representative of the actual acceptance for 
isotope separated fragments in the low energy runs. As can be seen the reconstruction in 
this case is close to perfect and purely statistical errors can be neglected. 



4-2 Tests of passage to 4tt using IQMD 



The following example shows the application to a case which cannot be a thermal distribu- 
tion since it is the superposition of the three isotopes (p, d, t) of hydrogen. As can be seen 
from Fig. [21 there are three separate contributions to the phase space distribution after 
applying the filter: they are from the CDC, the Helitron and the Zero Degree detectors. 
We show both the experimental data and the simulation with IQMD. Pion contributions 
are excluded or subtracted using the information obtained in 



The reconstruction of the full acceptance rapidity distribution of the simulation causes no 
difficulty as demonstrated in Fig. [31 The existence of low pt data due to the Zero Degree 
proves to be a very useful constraint on the reconstruction. 

The acceptance for identified protons in the higher energy runs was more restricted since 
the Zero Degree detector did not allow isotope separation and the need to use the Helitron 
in addition to the Plastic Wall restricted the forward part of the detector somewhat. This 
is shown in Fig. [H Still, we find that the reconstruction, for central collisions (which are 
of prime interest here) is very satisfactory. 

A case with a typical 'isotope acceptance' (for deuterons here) in the low energy runs, where 
only the CDC was available for isotope separation, is shown in Fig. [51 With the sharp filter 
distorted phase space cells near yo = ^1 due to multiple scattering effects in the target are 
also cut out. Due to missing low pt information and despite the mentioned constraints from 
the more complete data on Z = 1 fragments, the reconstruction has to allow for about 
7% errors indicated by the plotted systematic error bars. The accumulated counts in the 
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Fig. 1. Rapidity distributions of deuterons in Au+Au reactions at 0.25A GeV obtained from a 
thermal model with T = 77 MeV. The original predictions are given by (red) open squares, the 
reconstructed distributions after filtering and then extrapolating are given by (black) filled dots. 
The upper left panel shows the filtered distribution in scaled momentum versus rapidity space, 
the various color tones corresponding to yields differing by a factor 1.5. The lower panels show 
the integrated (left) and the constrained transverse rapidity distributions, while the longitudinal 
distribution is plotted in the upper right panel. 

experiment were actually 10 times higher than in the shown IQMD simulation, so that 
uncertainties originating here also from finite statistics were smaller in the experiment. 
On the other hand, again, the relatively faithful Air reconstruction is limited to central 
collisions, where low pt spectator material is not so copious. Here the sum-of-isotopes 
check proved to be mandatory to make sure there was no 'extrapolation catastrophe'. 
Constraints from charge separated data obtained with the ZERO Degree Detector, covering 
low transverse momenta, were used for these checks. 



4-3 Heavy clusters 



The 47r reconstruction of fragments with Z > 2, that we call here 'heavy clusters' (he) is 
somewhat problematic due to the more limited acceptance, although the two-dimensional 
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Fig. 2. Distributions dN/dutodyo of hydrogen fragments emitted in central collisions of Au+Au 
at 1.5A GeV. The various color tones correspond to cuts differing by factors 1.5. The distributions 
are filtered by a sharp-cut filter the left panel showing the experimental data, while the right 
panel shows the simulation (applying the same centrality and filter cuts) with IQMD. 

extrapolation procedure is an improvement over the conventional one-dimensional extrap- 
olation with separated rapidity bins. The results presented here for he, see also ref. [31], 
should be considered to be estimates when applied for example to 90° spectra. In Fig. |6] 
we show an illustrative picture of what it meant to 'reconstruct to 47r'. A stringent test 
using IQMD was not possible here due to limited statistics available for he in the simu- 
lation. However, observables, such as stopping, relying on in coverage, were found to be 
in very reasonable agreement with INDRA data [33j . As can be seen by inspecting Fig. [6l 
the increased 'stopping' at an incident energy of 0.25A GeV as compared to O.OQA GeV 
(larger transverse momenta relative to the longitudinal ones) is already suggested before 
the extrapolation. 
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Fig. 3. IQMD simulations: Longitudinal (upper panel), transverse and constrained transverse 
rapidity distributions of hydrogen fragments emitted in central collisions of Au+Au at 1.5 A GeV. 
The original predictions are given by (red) open squares, the reconstructed distributions after 
filtering (see Fig. [2]) and then extrapolating are given by (black) filled dots. 

4-4 Adjusting the low-energy runs to high energy runs 



As mentioned in section [21 the CDC was operated in split mode in the low energy runs, 
i.e. with a different voltage, [24], on the inner wires. Despite careful calibration to take this 
into account, it turned out by comparison with later runs made without the split voltage 
feature, that an additional correction was needed for the low energy runs. Fortunately this 
correction could be shown to involve just a simple rescaling factor. This is illustrated in 
Figs. [7] and [8] which consist both of two panels, the left one showing for central Au+Au 
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Fig. 4. Rapidity distributions of protons in Au+Au reactions at \.hA GeV obtained from 
IQMD-SM. The original predictions are given by (red) open squares, the reconstructed dis- 
tributions after filtering and then extrapolating are given by (black) filled dots. The upper left 
panel shows the filtered distribution in scaled momentum versus rapidity space. 

collisions at 0.4y4 GeV differences of the two operating modes (the earlier one with split 
voltage) before, and, in the right panel, after correction with a constant rescaling factor. 
Fig. [7] shows kinetic energy spectra for protons, while Fig. [8] shows transverse rapidity 
spectra for tritons. A common sharp filter is applied to the data from both experiments. 
The two peak-structure in Fig. [H] is a consequence of the limited acceptance, see the 
example in Fig. O The rescalings in both figures [7] and [H] was primarily along the abscissa, 
the rescaling of the ordinates then followed from the condition that the integrated yields 
were not affected. 

To avoid too many arbitrary parameters, the rescaling factors were fixed for all the low 
energy data down to 0.09A GeV. An example of the resulting extremely smooth behaviour 
of excitation functions for average kinetic energies covering the complete data is shown 
in Fig. O As will be shown later, this surprisingly regular trend is also predicted by 
IQMD simulations. We use scaled energy units here: a value equal to one indicates an 
energy equal to the incoming cm. kinetic energy per nucleon. The figure also shows the 
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Fig. 5. Rapidity distributions of deuterons in Au+Au reactions at 0.15 A GeV obtained from 
IQMD-SM. The original predictions are given by (red) open squares, the reconstructed distribu- 
tions after filtering and then extrapolating are given by (black) filled dots. The upper left panel 
shows the filtered distribution in scaled momentum versus rapidity space. 

reasonable agreement with data from ref. [12] which were obtained by our Collaboration 
using a different method. While we show here only a comparison for protons for reasons 
of space economy, similar conclusions hold for all other isotopes of hydrogen and helium. 



4 ■ 5 Z-renormalization 



By Z-renormalization we address the sum-of-isotopes check and correction after Att recon- 
struction, which uses the fact, shown earlier, that we had a more complete acceptance 
for fragments identified by charge only. This was most important for the low energy runs. 
We illustrate the necessary 'renormalization' in Fig. [10] for hydrogen and helium isotopes 
in central collisions of Au+Au at 0.15 A GeV. While renormalization is seen to be only a 
small effect for hydrogen, the worst case scenario shown for helium requires a significant 
correction for low transverse velocities, which can be traced back to errors due to not 
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Fig. 6. Example of 47r reconstruction of tlie phase space population for Z = 3 fragments in central 
(6o < 0.15) Au+Au collisions at two indicated incident energies. The upper two panels show the 
measured part after applying some smoothing, a sharp filter cut and reflection symmetry, the 
lower two panels show the result of a two-dimensional extrapolation to 47r geometry. 



Au+Au 0.4A GeV bo<0. 15 



lO'' 



CO 
Q. 
T3 



10^ 




Fig. 7. Left panel: kinetic energies (scaled units) of protons emitted at 90° in central collisions of 
Au+Au. Comparison of two experiments performed in the low energy (full black circles) and the 
high energy campaigns (red squares). In the right panel the low energy experiment is rescaled 
along the abscissa by a factor 1.1 to achieve agreement with the high energy experiment. 



accounting properly for Coulomb repulsion effects on low momenta not covered by our 
data. 
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Fig. 8. Sharp filtered transverse rapidity distributions for tritons in central Au+Au collisions 
at £^7^4 = OAA GeV. Comparison of two experiments performed in the high energy (open red 
squares) and the low energy campaigns (full black circles). Left panel: original data. Right panel 
with the low energy data rescaled by factor 1.1 along the abscissa. 
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Fig. 9. Average kinetic energy (scaled units) of 
protons emitted around 90° (cm.) in central col- 
lisions of Au+Au as a function of incident beam 
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5 Rapidity distributions 



We present the phase space distributions in a way which deviates somewhat from the tra- 
ditional way: instead of showing standard (longitudinal) rapidity distributions on a linear 
ordinate scale and then switching for the transverse direction to transverse momentum 
distributions on a logarithmic scale, we show the transverse degree of freedom in the same 
representation as the longitudinal degree of freedom, namely in terms of transverse rapidi- 
ties plotted on a linear scale as well (see definitions in section [3]). This allows to compare 
more directly longitudinal and transverse directions, a recurrent theme in the present 
work which is much concerned with assessing the degree of stopping and equilibration in 
these complex central collisions. We shall also be interested in the constrained transverse 
rapidity distributions (defined in section |3]) which are likely to be closest to thermalized 
distributions. 

An phase space distributions of five 'light charged particles' (LCP: p, d, t, ^He,''He) have 
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Fig. 10. Left panel: Constrained transverse rapidity distribution of H fragments in central Au+Au 
collisions. (Black) open squares: demanding Z identification only, (pink) closed circles: sum of 
identified isotopes (p, d, t). Right panel: same for He. 

been reconstructed using four different centralities in 25 system energies, i.e. 500 2d- 
spectra are available. This number could be doubled if we consider that two different 
selection criteria (ERAT and MUL, see section [3]) were applied throughout. However, 
unless otherwise stated, we shall show generally only the ERAT selected data. Further, for 
E/A < 0.4 GeV we have the data for heavier charge separated clusters (typically up to 
Z=8) and for E/A > 0.4 GeV charged pion data of both polarities (see [SHj). It is out of 
question to present here all this rich information. Instead we shall present samples which 
serve to illustrate some of the typical aspects of the LCP data and add a few remarks 
on heavier clusters and pion data [38]. In later sections some of these aspects will be 
summarized in terms of simple concepts: 'stopping', 'radial flow' and 'chemistry'. These 
concepts reflect the time sequence of the central reactions. Earlier work on these aspects 
will be referred to. 

Some of the rapidity distributions will be compared with distributions expected from a 
Boltzmann thermal model. These thermal distributions are generated varying T until the 
variance of the constrained transverse rapidity distribution (i.e. with the scaled longi- 
tudinal distribution constrained by \yzo\ < 0.5 and within the yxmo range shown in the 
respective panels) is reproduced (the same constraints are used on the thermal simulation). 
The obtained T will be dubbed 'equivalent' temperatures, Teq. The longitudinal constraint 
leads to a cutoff of part of the yields, a cutoff which is particle-mass and T dependent. 
This cut is corrected for when extracting so called effective 'mid-rapidity yields' from the 
experimental data by comparison to the thermal model (see section [8]). 



5.1 Light charged particles (LCP) 



The next two 4-panel figures. Figs. [HI [T21 show the three kinds of rapidity distributions 
for protons emitted in the most central Au-(-Au collisions at two incident energies differing 
by one order of magnitude: 0.15^4 GeV and 1.5^4 GeV. 
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Fig. 11. Upper left panel: Longitudinal (blue open triangles) and transverse rapidity distribu- 
tions (red full diamonds) of protons in central collisions of Au+Au at 0.15 A GeV beam energy. 
The third, lower yield curve (magenta open diamonds) in the panel represents the 'constrained' 
transverse rapidity distribution: here a cut \yzo\ < 0.5 on the longitudinal rapidity is applied. 
The smooth curves are just guides for the eye here. In the other three panels the three kinds of 
distributions are compared to a thermal distribution (corresponding to an equivalent tempera- 
ture of 43 MeV) having the same area as the respective data, but a variance as implied by the 
constrained transverse distribution (seen in upper right panel). 

The upper-left panels in the two figures show that the longitudinal rapidity distributions 
are broader than the transverse distributions, but while the effect is relatively modest at 
the lower energy, it is more conspicuous at the higher energy. The constrained transverse 
distributions are generally somewhat broader than the integrated transverse distributions 
and can be reproduced rather well by a thermal model simulation having an adjusted 
'equivalent' temperature, T^g, and the same total area. We find T^g = 43 MeV, and 158 
MeV, respectively, for protons at the two energies. This Teg comes closest to the 'inverse 
slopes' usually fitted to pt spectra in the literature, but is different in the sense that it is 
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Fig. 12. Upper left panel: Longitudinal (blue open triangles) and transverse rapidity distribu- 
tions (red full diamonds) of protons in central collisions of Au+Au at 1.5^ GeV beam energy. 
The 'constrained' transverse rapidity distribution is also shown: here a cut \yzo\ < 0.5 on the 
longitudinal rapidity is applied. The data points are joined by smooth curves to guide the eye. In 
the other three panels the three kinds of experimental distributions are compared to a thermal 
distribution (smooth black curves) having the same area as the respective data and a variance 
as implied by the constrained transverse distribution (seen in upper right panel). 



obtained by just demanding a reproduction of the bulk experimental variances taken in 
the shown rapidity ranges and hence excluding far out tails of the distributions. On the 
linear scales applied in the figures one finds a very satisfactory reproduction of the shape of 
the dN/dUxmo distributions (upper right panels). A comparison to the equivalent thermal 
representations fixed from the constrained data, to the full transverse and longitudinal 
distributions in the lower panels stresses the differences to the naive thermal expecta- 
tions. The full transverse distributions are modestly narrower only, but the longitudinal 
distributions clearly deviate. 
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Fig. 13. Central (6o < 0.15) Au+Au collisions at 0.4A GeV. Comparison of experimental scaled 
and constrained transverse rapidity distributions dN/dyxmO of protons (upper right), deuterons 
(lower left) and tritons (lower right) with thermal distributions (smooth black curves) assuming 
the indicated equivalent temperatures T^q. The upper left panel shows the distributions of the 
three H isotopes, their sum (pink full dots) and the distribution for Z=l (excluding created 
particles) obtained requiring only charge identification. 

The next two 4-panel figures, Figs. [131 El take a closer look at constrained transverse 
rapidities varying the (hydrogen) isotope mass. The system is Au+Au at OAA GeV with 
bo < 0.15. First we show the data, then a simulation with IQMD. In each case we make 
in the upper left panel the sum-of-isotopes check (which is fulfilled by definition in the 
simulation) and we compare in the other panels with a thermal model calculation adjusting 
Teg. These one-shape parameter fits are close to perfect, but require a rising equivalent 
temperature with the isotope mass, a well known phenomenon often interpreted as 'radial 
flow' in the literature that we shall come back to in section [71 Clearly there is no global 
equilibrium, but at very best a 'local' equilibrium (in the hydrodynamic sense). 
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Fig. 14. Central (6o < 0.15) Au+Au collisions at 0.4A GeV. Comparison of simulated 
(IQMD-SM) scaled and constrained transverse rapidity distributions dN/dyxmO of protons (up- 
per right), deuterons (lower left) and tritons (lower right) with thermal distributions (smooth 
black curves) assuming the indicated equivalent temperatures T^q. 

A closer look at the simulated data reveals small but systematic shape differences that 
exceed those found when using experimental data instead of the IQMD events. One could 
argue that the experimental data look more 'thermal'. Note the rather strong decrease of 
yields with isotope mass in contrast to the more comparable yields in the experiment shown 
in the previous figure. Also, if one were to associate, naively, the strength of radial fiow with 
the mass dependence of T^g, then one would conclude that the simulation overestimates 
the fiow somewhat, since, as indicated in the various panels the Tgg rise from 85 MeV for 
protons to 116 MeV for tritons in the data, while for the simulation we find 67 MeV and 
118 MeV, respectively. For further discussions of the T^g values, see later section [71 

The Fig. [15] showing constrained transverse rapidity distributions for central collisions in 



21 



Ru+Ru 1 .5A GeV bo<0.1 5 He 




-1.0 -0.5 0.0 0.5 1.0 



12 



T3 



deuteron 



T =179MeV 

eq 




L-..-al" I I I I I I I "IB--. 

-1.0 -0.5 0.0 0.5 1.0 




T3 



0.8 



0.6 



0.4 



0.2 



I ' r 



yxmO 

' " T„ =234 MeV 

/ eq 

p-like 




-0.8 -0.4 0.0 0.4 0.8 

yxmO 



-0.8 -0.4 0.0 0.4 0.8 

yxmO 



Fig. 15. Comparison of experimental constrained transverse rapidity distributions dN/dyxmO 
with thermal distributions (black smooth curves) having the same first and second moments in 
central collisions of Ru+Ru at l.hA GeV. The various ejectiles and the equivalent temperatures 
T(,q are indicated. In the lower right panel we also plotted the ^He distribution (dashed smooth 
curve, 'p-like') expected if the equivalent temperature was the same as for protons, but with the 
integrated yields unchanged. 



the system Ru+Ru at 1.5 A GeV serves two purposes: 1) besides hydrogen isotopes (p and 
d) it also shows data for ^He and '^He in the lower panels. Again thermal model calculations 
with the indicated T^q are shown along with the data; 2) We show that the ^He data are 
not well represented by a coalescence model using the measured proton data as generating 
spectrum: see the dashed curve. Indeed the observed rising equivalent temperatures with 
mass contradict the naive model. This is at variance with ref. ^3] where the validity of 
the power law for a wide range of incident energies and centralities in Au+Au systems 
was stressed. The authors mentioned that a cutoff of pt/-^ smaller than 0.2 GeV/c was 
required and that the single ('free') proton spectra had to be used, rather than the total 
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Fig. 16. Central (6o < 0.15) (bottom) and semi-central (top) Au+Au collisions at 0.4A GeV. 
Comparison of experimental scaled and constrained transverse rapidity distributions dN/dyxmO of 
^He (left) and ^He (right) fragments with thermal distributions (smooth black curves) assuming 
the indicated equivalent temperatures T^q- 

spectra including those bound in clusters as one would naively expect. Non-perturbative 
features of clusterization are suggested strongly by our data, as we shall show later in this 
section and in section [8l 

In the cases shown so far we observe a regular rise of Teq with the mass of the ejectile. 
However, frequently in the literature, the so-called 'helium anomaly' is mentioned |44j . 
namely the observation that ^He kinetic energies are not lower, but higher than the "^He 
energies, a phenomenon that to our knowledge has not been reproduced by microscopic 
dynamic reaction models. In Fig. [11] we show for Au+Au that even at the relatively high 
incident energy of E/A = 0.4 GeV this phenomenon still subsists to some degree as the 
Teq for ^He are not found to exceed those of ^He. Further we show in the panels the rather 
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Fig. 17. Left panel: Constrained transverse rapidity distributions of protons in central collisions 
of Au+Au. Comparison of data (red closed diamonds) and IQMD simulation (blue open circles). 
Right panel: comparison of experimental data for ^He (red closed diamonds) and ^H (blue open 
diamonds) in the same reaction. 

strong drop of Teq with decreasing centrality, which is varied from bo < 0.15 (lower panels) 
to 0.45 < &o < 0.55 (upper panels). This is a general feature also observed at other energies 
and was reported earlier 



5.2 The influence of clusterization 



The failure of the naive coalescence approach consisting in trying to reproduce the spectra 
of heavier clusters from the measured nucleon (proton) spectra (Fig. [T5l) and some of the 
'He-anomaly' discussed before may be connected with our current failure to understand 
clusterization quantitatively on a microscopic level (compare also the upper right panels 
of Figs. [T3] and [T^ . Concerning the probability of nucleons to cluster, we are at SIS in a 
non-perturbative regime (see section [8]): heavier cluster formation has a back-influence on 
the lighter generating transverse rapidity spectra. 

The two-panel Fig. [17] illustrates the non-perturbative features at SIS. In the left panel 
we show the remarkable difference between FOPI and IQMD for the transverse rapidity 
distributions of protons. The surplus of IQMD for lower transverse velocities or momenta 
is due to a lack of sufficient clusterization: in the experiment more copious cluster for- 
mation massively depletes the low momenta. The right panel compares two experimental 
distributions: now it is the ^He that appears to have its low momenta depleted relative 
to ^H (the naive perturbative coalescence model does not predict any difference here). In 
view of the finding of the left panel, it is tempting to associate the effect to the formation 
of heavier clusters from the 'primeval' ^He (created in earlier expansion stages). 

This conjecture is supported by Fig. [18] which shows the '^H-'^He difference spectrum to- 
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gether with data for ^He. The ^H and the "^He compete to be a condensation nucleus to a 
possible ^He. If both mass 3 isotopes are in a neutron- rich environment, the ^He will 'win' 
for two reasons: 

a) it is easier to 'find' a single neutron to attach to ^He than a single proton to attach to 
3H; 

b) in contrast to ^H, the ^He nucleus does not Coulomb-repulse its needed partner. 

A quantitative transport model theory must include the formation of a clusters if these 
conjectures are correct. 



5.3 How to define stopping: varxz 



The final figure in this section on rapidity distributions, Fig. [T^l serves to illustrate what 
we shall call 'partial transparency' or 'incomplete stopping'. Transparency increases as 
either the studied system mass is lowered, or the incident energy is raised by going from 
E/A = 0.4 GeV (lower panels) to 1.5 GeV (upper panels). With this figure, which shows 
data for the most central collisions {bo < 0.15 or 2.2% of the equivalent sharp cross 
section), we argue that a good measure of the degree of stopping consists in evaluating 
the ratio varxz of variances varx and varz of the transverse relative to the longitudinal 
rapidity distribution, respectively, both of which are shown in the figure for deuterons. 
Note also the subtle evolution of shapes, especially at the higher energy: there is no trivial 
subdivision into 'participants' or 'spectators' possible. The interpretation in terms of two 
remnant counterflowing (not completely stopped) 'fluids' is strongly supported by isospin 
tracer methods (see section |6]) and suggests that not only global, but also local equilibrium 
is not achieved. 
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Fig. 19. Comparison of scaled rapidity distributions for deuterons in central collisions at 
EjA = 1.5 GeV (top panels) and 0.4 GeV (bottom panels): shown are data for Ca+Ca, Ru+Ru 
and Au+Au (from left to right). In each panel, the scaled longitudinal rapidity y^o distributions 
(blue open triangles) are shown together with the scaled transverse rapidity y^o distributions 
(red). The derived value of the stopping observable varxz is also indicated. 

6 Stopping 



As introduced in section [5l3| we quantify the degree of stopping by comparing the variances 
of transverse rapidity distributions (defined in section [3]) with that in the longitudinal (2;) 
direction. Throughout we use scaled rapidities (i.e. the rapidity gap is projected unto the 
fixed interval (—1, +1)) and typical samples have been shown and discussed already (sec- 
tion E]). All rapidities are evaluated in the cm.. However, for the ratios we discuss the 
scaling is immaterial. We shall call varxz the ratio of the variances in the x-direction and 
the z direction. The choice of x is arbitrary, i.e. it is not connected with the azimuth of the 
reaction plane, varxz(l) is defined like varxz except that the integrations for calculating 
the variance are limited to the range \yzo\ < 1; in [32j this was called vartl and \yxo\ was 
also constrained to be < 1. Such a restricted measure of stopping can be more adequate 
if measured data outside the range \yzo\ < 1 are missing or less reliable. For varxz(0.5) 
the transverse rapidity distribution is taken under the constraint \yzo\ < 0.5 on the scaled 
longitudinal rapidity. This measure of stopping is clearly biased and has systematically 
higher values, but is of interest when comparing with many data in the (higher energy) 
literature, where often attention is focused on transverse momentum spectra around (lon- 
gitudinal) midrapidity. When no number in parenthesis is given, the intervals in rapidity 
are sufficiently broad to make the value of varxz asymptotically stable within the indicated, 
mostly systematic, errors. 
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There are other ways of characterizing stopping. Videbaek and Hansen [16] have introduced 
the mean rapidity shift < 6y > which is defined by 

.. J'^Z\y-ytgt\{dNp/dy)dy 

y - P_Z{dN,/dy)dy 

where ycm is the rapidity of the cm., y^gt is the target rapidity and dNp/dy is the proton 
rapidity distribution (a similar formula exists relative to the projectile rapidity yp with 
the integration limits interchanged). This observable was also determined in some works 
of our Collaboration, [ISPU] . The inelasticity K was used by NA49 [17] and is defined 
by i^ = Einei/iy/s/l — nip) with Emei being equal to the average energy loss of incident 
nucleons with rest mass nip. The common feature of these definitions is that one assesses 
stopping relative to initial conditions: the conclusion is then that a high degree of stopping 
is reached in colhsions all the way up to at least SPS energies, especially if it is assessed 
in terms of energies. Our stopping variable has a different goal: we want to assess the 
difference relative to a completely stopped scenario, an assumption that is frequently 
assumed in (one-fluid) hydrodynamic codes for heavy ion collisions because of its inherent 
simplicity of a needed starting point of the calculation. Our observable, varxz, is important 
to assess the partial memory of the original (accelerator induced) counterflow of two fluids 
and hence deviations from local equilibrium. It is also of value to help determine viscosity 
properties of nuclear matter [50] • 



6.1 Isospin tracing 



We invariably flnd that longitudinal rapidity distributions are broader than transverse 
rapidity distributions, see section O In principle, for a given system, one cannot exclude 
that this phenomenon is caused by a rebound opposite to the incident direction after 
a complete stop. The system size dependence of the effect, see namely Fig. [191 can be 
used to strongly argue against this interpretation which would imply an unlikely stronger 
rebound in smaller systems. An alternative method to demonstrate partial transparency 
was introduced by our Collaboration [52] and consists in so called isospin tracing. We refer 
the reader to the original publication for further explanations. Briefly, we combine rapidity 
distribution information of four systems involving the isotopes ^^Zr and ^^Ru : Ru+Ru, 
Zr+Zr, Ru+Zr and Zr+Ru, where the mixed systems are merely technically different in 
the sense that beam and target are inverted. 

In the following, abbreviating N^'^^^{p) the rapidity distribution of protons in the reaction 

Ru+Zr, etc, we deflne 

R^y(p) = [N^^^'-ip] - Ar^^^^"(p)]/[iVy^«^"(p) - N^'-^'^ip)] 

This observable is a slight variant of the observable Rz used in [52]: it is an average of 

the two opposite-sign branches obtained by switching target and projectile in the mixed 

system. An arbitrary (here negative) equal sign is deflned (the two branches of Rz in the 

chosen symmetric- mass systems must be equal except for sign). In Fig. [20] we show -R4y 

for central collisions at 0.4A and 1.5^4 GeV together with rapidity distribution plots. In 

case of complete mixing with loss of memory (both alternatives: rebound as well as partial 



27 



T3 



2*1 


1 ' 1 








Ru+Ru 


A 


♦ transverse . 




0.4A GeV 


A longitudinal 


20 


"b0<0.15 


n 


- 


16 


- 


i 


I 






^ 


I. 


12 


- 


.5 


i. 






A? 


^ 


8 


t 


? 


^A 




* 




* 




$ 




$ 


4 


♦ 




♦ 




ft 




A 




, ./ , 




. , ^, , " 





Ru+Ru 




{} 


' 1 ' 

♦ transverse 


60 


-I.SAGeV 
- b0<0.15 




A 


A longitudinal 


50 


— 




1 I 


— 


^40 


- 




fA'^^^^AJ 


* — 


i 30 




A^ 




^A : 






'^J 




s^ 




A 






A 


20 


A 


♦ 




♦ — 
A 

♦a 


10 


A ♦ 






♦ A 

, r»^A , 



cc 




■si- 
ce 




Fig. 20. Top panels: Longitudinal and transverse rapidity distributions of protons in collisions 
of Ru+Ru at OAA GeV (left) and 1.5A GeV. Bottom panels: The corresponding R^y observables 
for two slightly different indicated centralities as a function of the longitudinal rapidity yzo- 

transparency) of the incident geometry, the rapidity dependences of R4y should be fiat at 
zero value. There is indeed a small flat part around i/zo = 0, especially for the most central 
{bo < 0.15) selection, but for higher \yzo\ there is an increasing deviation from zero, the 
sign of which can be unambiguously associated with transparency, rather than rebound. 
The effect is definitely more pronounced at the higher energy in full accordance with the 
information from the variances of the rapidity distributions (see the upper panels in the 
figure). Only statistical errors are shown in the lower panels of Fig. [201 A more detailed 
presentation and discussion of the new isospin tracing data will be published elsewhere 

ESI. 



6.2 Global stopping and the EOS 



By adding up the measured rapidity distributions of the different ejectiles one can try to 
obtain a global system stopping. The excitation function for Au on Au, first published 
in [32], was later slightly revised and combined with data from the INDRA collaboration 
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Fig. 21. Excitation function of tlie stopping observable varxz(l) in central Au+Au collisions 
(6o < 0.15) |32|33j . FOPI data: open blue triangles, INDRA data: black dots. These experimental 
data points are joined by straight line segments to guide the eye. IQMD simulations extending 
from 0.12A to 1.5^ GeV: red full triangles (HM) and green full diamonds (SM). 



|33j . The combined data are shown again in Fig. |2TJ but this time together with IQMD 
simulations using alternatively a hard {K = 380 MeV, HM) or a soft {K = 200 MeV, 
SM) EOS with momentum dependence [^. As can be seen from the Figure, the difference 
between the HM and the SM simulation is significantly larger than the error bars in 
the data. This is an important observation as it opens up the possibility to observe EOS 
sensitivity also in central collisions in contrast to directed and elliptic flow which are linked 
to off-centrality (and in general collisions with less achieved compression). However, it is 
also clear that the 'residual' interaction, i.e. the explicit collision term, influences the 
outcome. The present parameterization of IQMD as used here is obviously not able to 
reproduce the data, in particular the rapid drop of varxz(l) beyond 0.8A GeV is not 
reproduced. A fair reproduction of a portion (0.25A to l.OA GeV) of the excitation function 
was achieved in 



Clearly, varxz data should prove useful also to constrain the viscosity of nucleonic (and 
more generally of hadronic) matter if compared to transport model simulations. It is 
tempting to conjecture that the rising part of the excitation function is dominated by 
the decreasing importance of Pauli blocking as the rapidity gap exceeds the Fermi energy. 
Particle creation starts becoming important beyond lA GeV [38] (see also our Fig. 41 in the 
present work). But the decreasing trend of the varxz data seen in Fig. HH suggests that this 
new channel does not compensate for the decrease and forward focusing of elastic nucleon- 
nucleon processes. The rapid decrease of varxz at the high SIS energy end is continued at 
AGS and SPS as we shall show later for net protons. 
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Fig. 22. Excitation function of the stopping observable varxz in central Au+Au collisions 
{bo < 0.15). As indicated in the various panels the excitation functions are shown separately 
for the three hydrogen isotopes. The right-lower panel allows to compare relative to the average 
weighted Z=l fragments. The smooth curve is just a polynomial fit. 

6.3 p, d, t stopping: a hierarchy 



By plotting varxz separately for identified fragments further insights (and constraints to 
simulations) can be gained. This is shown in Fig. [221 A remarkable feature is the similarity 
in the behaviour for the three isotopes, that suggests the collectivity of the phenomenon. 
There is a well defined maximum around 0.4^4 GeV for each of the isotopes qualitatively 
similar to the observed [32] global stopping. On the other hand a closer look reveals that 
there are some subtle differences in the shapes of the excitation functions. Except around 
the maximum, there is a hierarchy in the degree of stopping. This is illustrated in Fig. | 



On the average clusters have 'seen' less violent (less stopped) collisions when the stopping 
is incomplete for single emitted nucleons (away from the maximum around 0.4^4 GeV). 
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Fig. 23. Mass number dependence of the stopping observable varxz of hydrogen fragments in 
central Au+Au collisions {bo < 0.15). Data are shown for 4 representative energies indicated in 
the panels (notice the one order of magnitude span) . For comparison data are also shown for the 
'constrained' stopping observable varxz(0.5) (blue open triangles, error bars omitted). The lines 
are linear fits. 



This effect may be a surface effect but is not a 'trivial' spectator-matter effect, as the 
rapidity distributions in these most central collisions do not exhibit well defined 'spectator 
ears' around \yzo\ = 1: there are no spectators in the strict sense here. The constrained 
stopping varxz(0.5) (blue open triangles in the figure) is always somewhat larger because 
the transverse rapidity distributions are broadened when a cut around longitudinal rapidity 
is applied. The trends with mass number and incident energy are similar however. 
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6-4 System size dependence of stopping 



As remarked earlier, the system size dependence of the observable varxz is an important 
constraint to the question of transparency versus rebound dominated collision scenario, es- 
pecially when only the most central collisions are compared in scaled units of centrality. All 
our data support transparency dominance. This does not mean that the degree of stopping 
is small, but it is significantly less than expected from ideal one-fluid hydrodynamics. 

Figure [2l] shows a rather complex behaviour and will be a challenge to ambitious mi- 
croscopic simulations. One remarkable feature is the apparent saturation of varxz at a 
relatively low value for the highest incident energy. This was already seen earlier, |32] , for 
global stopping. This suggests the beginning of a phenomenon seen in a more spectacular 
way recently at the SPS [47j {E/A = 158 A GeV): at midrapidity there is no non-trivial 
increase of population in Pb-|-Pb relative to p-|-p collisions. 



6.5 Isospin and stopping 



One expects isospin dependences because the free nucleon-nucleon cross sections ann and 
anp are different and also because the mean fields, shown earlier to influence global stop- 
ping, are expected to be isospin dependent. However the data show no convincing effect 
considering the error limits. This is shown in Fig. [25] where a comparison of two sys- 
tems with the same mass, but different composition (Ru+Ru, N/Z = 1.182 and Zr+Zr, 
N/Z = 1.400) is shown. In this context it is interesting to note that very recently a similar 
observation was made [5Tj at lower energies (32 — lOOA MeV) for Xe-|-Sn using various 
isotopes. 



6. 6 Centrality dependence of stopping 



For comparisons with simulations it is useful to check how sensitively the observable varxz 
reacts to matching the centrality correctly with that of the experiment. The relatively 
strong dependence of varxz on centrality b^ is illustrated in Fig. [22] for Au+Au at two 
rather different incident energies. The linear fits allow to estimate a limit for 6o = 0. Since 
varxz is the ratio of varx and varz (the zero suffix is immaterial here since the scaling 
cancels out ), it is interesting to explore what causes the strong centrality dependence. 
The figure also illustrates the comparatively modest change of the scaled variance, varxmO 
of the constrained transverse rapidity distribution (which we remind is obtained with a 
with a cut I2/20I < 0.5 on the scaled longitudinal rapidity). 



6. 1 Heavy clusters 



Since at lower energies we have been able to get information also for fragments with Z > 2, 
we can extend the hierarchy of stopping to heavy clusters. The results are put together 
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Fig. 24. System size dependence of varxz for protons (black circles), deuterons (red squares) and 
tritons (blue triangles) at various indicated beam energies. 



for three incident energies and three systems in Fig. [271 This confirms the existence of a 
stopping hierarchy in a very systematic way. Also a significant system-size dependence is 
evidenced. Note that the Z-dependence is flattest for Au+Au at 0.4^4 GeV, close to the 
maximum of global stopping shown before in Fig. |2TJ An extension, for O.ISA GeV, all 
the way to Z=20 using INDRA data has been shown earlier in |33], demonstrating also a 
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Fig. 25. Mass dependence of varxz for hydrogen isotopes in central (6o < 0.15) Ru+Ru and 
Zr+Zr collisions at incident beam energies of OAA GeV (left panel) and 1.5^ GeV (right panel). 
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Fig. 26. Centrality dependence of proton stopping (varxz) and mid-rapidity constrained scaled 
transverse variance (varxmO) at 0.15 and 1.5A GeV. 

reasonable consistency between data of FOPI and INDRA. Quantitatively, such data are 
presently outside the capabilities of microscopic simulations, but present a challenge for 
our theoretical understanding. 

6.8 Joining up to higher energies 



As mentioned at the beginning of this section, stopping at higher incident energies has 
not been characterized by varxz. The reason is connected with the increasing difficulty 
to measure and reconstruct the distributions dN/dyz and dN/dyx over the full range of 
rapidities that one needs to know. Often transverse momentum spectra in narrow yz bins. 
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Fig. 27. Charge number and system size dependence of the stopping observable varxz for frag- 
ments emitted in central Au+Au collisions (6o < 0.15). Data are shown for 3 representative ener- 
gies and three systems indicated in the panels. The smooth curves, fitted to oq + a-iZ~^ + aiZ , 
serve to guide the eye. 

preferably around ly^ol = 0, are fitted individually to thermal ansatzes. We show therefore 
the ratio of variances using in the numerator the variance inferred from inverse slopes 
Tinv determined from transverse momentum spectra in a bin Ay^o around (longitudinal) 
mid-rapidity. In terms of the notation introduced at the beginning of this section, this is 
likely to give an upper limit, varxz{r^ 0.1), to the unconstrained varxz, as suggested by 
Fig. I 



With these restrictions in mind, we show in Fig. [28] an extension of our proton data on 
stopping to higher energies. The AGS (Alternating Gradient Synchrotron at Brookhaven) 
data point was constructed from ref. [54j|, the SPS (CERN Super Proton Synchrotron) 
data point from ref. ^55j (NA49). A rather dramatic decrease of stopping, as defined in the 
present work, is evidenced, suggesting a fair amount of residual counterfiow if one takes 
transverse fluctuations as scale. The straight line is a fit joining to SIS data starting at 
0.6A GeV. However for the interpretation of stopping excitation functions of one specific 
particle, here protons, extending over more than three orders in the magnitude of the 
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Fig. 28. Excitation function of the stopping ob- 
servable varxz for protons in central Au+Au col- 
lisions. In addition to FOPI data, obtained at 
SIS, we show here estimates from AGS (red full 
square) and SPS (blue full triangle) data. See 
text for further details. 
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incident energy, it is important to realize that they can be misleading on the global system 
trend. The reason is that stopping hierarchies are omnipresent. At the low energy end 
protons are the most stopped particles, while at the high energy end net protons (the 
difference between protons and antiprotons) tend to be rather the least stopped particles. 
We are lead to expect stopping hierarchies in the higher energy regimes similar to the ones 
described here, but shall not pursue this as it would far exceed the scope of the present 
work. As consequence of incomplete stopping, the time the system spends at maximum 
compression (and heat) is shortened. This could prevent the full development of critical 
fluctuations expected theoretically in certain areas of the phase diagram. 



7 Radial Flow 



One of the first attempts to understand radial flow in central heavy ion collisions in the 
framework of microscopic transport theory was published in [50] ■ There, the transverse 
flow energy El per nucleon was deflned by 



El = JJdrdpfilmivify JJdrdpf 



and the transverse flow velocity v{ was calculated locally, using only those particles from 
the environment that have participated in the collisions and stopping the integration when 
the density po/8 was reached in the expansion stage (/ is the Wigner distribution function). 
Evidently this kind of radial flow is not a direct observable. Rather, it is usually inferred 
from the mass dependence of average kinetic energies of particles emitted around mid- 
rapidity (in cylindrical coordinates) or around 90° cm. (in spherical coordinates). This 
implies an interpretation of the data suggested by hydrodynamical thinking: the system's 
quasi-adiabatic expansion leads to cooling, part of the energy being converted to (radial) 
flow. As we shall see this interpretation can be put to a critical test by studying the 
dependence of the apparent radial flow on the centrality and the system size and also by 
inspecting the correlation of the radial flow with the degree of clusterization taking the 
latter as a measure of the degree of cooling. 
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Fig. 29. Average kinetic energy per nucleon (scaled units) of protons, deuterons and ^He (left 
panel) and of ^H and ^He (right panel) emitted in 90° it 10° in the cm. in central collisions of 
Au+Au as function of beam energy. The excitation functions follow a linear dependence (straight 
lines) if plotted on a logarithmic energy scale. 

Relevant data for the SIS (or Bevalac) energy range have been published by Lisa et al. |15] 
and in a number of papers from our collaboration [^42f48f57f58] . Azimuthal dependences 
of radial flow in the same energy range were studied in refs. [281130] and shown to be 
sensitive to the EOS, while little sensitivity was suggested for the azimuthally averaged 
radial flow. 

We shall start by taking a look at the systematic evolution with incident energy and 
ejectile mass (hydrogen and helium isotopes) of average scaled kinetic energies per nucleon, 
Eko/A. This is shown in Fig. |29]for collisions with bo < 0.25 and emissions into cm. angles 
(90 ± 10°). A striking feature is the extremely smooth evolution with beam energy: the 
point to point errors seem to be smaller than the indicated systematic errors. In these 
scaled units the energies gradually decrease as the available energy must be shared by an 
increasing number of emitted particles. The regularity suggests a very continuous change 
without any dramatic fast transition to some higher entropy phase. Note that at the low 
end we observe protons emitted perpendicular to the beam axis with on the average twice 
the incident beam energy per nucleon. 

Another feature, demonstrated in the right panel of the figure is the difference between 
the two equal mass isotopes ^He and ^H. Starting at the low beam energy end with a 
difference of kinetic energies well in excess of the expectations from Coulomb effects [12] , 
this difference gradually converges to zero in these scaled units. We have already discussed 
this 'anomaly' in section [5] for the special case of an incident energy of 0.4^1 GeV and have 
given there a tentative interpretation. Here we see the 'anomaly' apparently disappearing 
as production of clusters heavier than mass three drastically decreases (see also Fig. H3|) . 



Another demonstration for the need to reproduce correctly the degree of clusterization in 
order to understand the kinetic energies of emitted particles is given in Fig. [30] for protons. 
Much like the experimental data the calculated data are reproduced with high accuracy 
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Fig. 30. Average kinetic energy (scaled units) of protons emitted around 90° cm. in central 
collisions of Au+Au as function of the beam energy. Comparison of experimental data with 
simulations using IQMD-SM. 

(0.7%) by a linear dependence on \og{E/A). The slope characterizing the calculated trend 
is different however. The reason is connected with a lack of sufficient clusterization at the 
low energy end in the simulation. Too many single protons, together with the necessity to 
conserve the total energy (and mass), lead to smaller energy per nucleon. This stresses the 
need to reproduce the degree of clusterization for quantitatively correct comparisons with 
experimental observables pertaining to specific identified particles. In particular radial fiow 
is connected with cooling leading to clusters. The two observables, cluster yields and radial 
flow, are therefore interrelated. 

In Fig. [31] we show average kinetic energies of H and He isotopes emitted at 90° cm. in 
Au on Au collisions with centrality bo < 0.25. The data, plotted versus fragment mass for 
various incident beam energies, show pronounced structures that tend to decrease with 
incident energy: by adding INDRA data [52] for Xe+Sn at 0.05 GeV a very large range 
of incident energies is spanned. The structures prevent a straight-forward determination 
of a linear slope with mass that could be tentatively associated to a collective radial flow. 
The straight lines plotted in the various panels join the proton and deuteron data to 
the average kinetic energy of ^He and '^H and ignore the ^He data. As can be seen, the 
average mass three value continues the proton-deuteron trend with a remarkable accuracy, 
an observation that we found to be true for all studied system-energies and centralities 
(see also Fig.EZj), allowing therefore to postulate a well deflned observable, the slope of the 
corresponding straight line, that we shall call 'radial flow with the A3 method', although 
this interpretation is clearly subject to caution. 

Deflning Ek = Ef x A + Eo with Ef = (7/ — 1) rriN {ttlm nucleon rest mass), we convert the 
mass slope Ej from the just introduced A3 method to a (radial) flow velocity /3j , derived 
from the average flow factor 7^ , and the offsets Eo in Fig. [HI] to an 'offset temperature' 
To using relativistic formulae and ignoring mean fleld contributions to E^, , in particular 
Coulomb flelds. This allows us to compare to the analysis of |15], although the authors 
seem to have handled the data structures differently. 
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Fig. 31. Average kinetic energies E^ of H and He isotopes emitted at 90° cm. in Au on Au 
collisions with centrality bo < 0.25. Full (blue) diamonds: ^'^He. The data are obtained at various 
indicated beam energies. The straight line is a linear fit up to mass 3. The data point at mass 
3 (x symbol) was taken to be the average of ^H and ^He. The derived apparent temperatures, 
To, and (radial) flow velocities, /?/, are indicated. The upper left panel illustrates the use of the 
same method to INDRA data [59] for Xe+Sn taken at 50^1 MeV. 



This is shown in Fig.|32l The two datasets with nearly matched centralities are compatible, 
although our offset temperatures are systematically somewhat higher. Our data extend 
over a larger energy range and follow a more regular trend fixed by smaller error bars. The 
interpretation of the data in terms of a collective flow and a common 'local' temperature 
is probably too naive: the stopping hierarchies evidenced in section [5] make a perfectly 
common flow improbable. We expect further the apparent offset To parameters to be 
influenced both by Coulomb fields (especially at the low energy end) and by the repulsive 
nuclear fields built up by the compression (see Fig. [2T]) . 



7.1 System size dependences 



In Fig. [33]we show the system size dependence of radial flow varying the incident energy by 
one order of magnitude. The increasing trend with system size is striking and characteristic 
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Fig. 32. Radial flow velocities /3/ (left panel) and offset temperatures To (right panel) for Au+Au 
collisions deduced from 90° kinetic energies of mass 1-3 ejectiles {A3 method). The present data 
(blue full symbols) are compared to data from ref. [45]. The smooth curves are fitted to the 
present data only and serve to guide the eye. 

for all energies. In each panel two sets of data are shown: the data points plotted with the 
open symbols and the larger error bars were obtained fitting two parameters separately 
to the data for one incident energy, the flow (slope) and the offset To. The second set of 
data points were obtained using a common offset To indicated in the panel, a constraint 
allowing to fix the flow value better. 

The difference between the two curves is small. However we find, when allowing two 
adjustable parameters for each system, that the To values is not completely independent 
of the system size: from the lightest system to the heaviest (Au+Au) there is an increase by 
about 20%. This could be a 'relic' of the repulsive mean fields due to increased compression 
in heavier systems. Unfortunately the uncertainty of To is relatively large. 



7.2 Comparison to IQMD 



The observed system size dependence of the radial flow is not trivial: a comparison with 
the simulation using IQMD-SM, see Fig.|3ll reveals that the strong dependence of the data 
on system size is not correctly reproduced. In |60], using IQMD, the authors conclude that 
radial flow is not sensitive to the EOS: they find that the sum of two contributions to this 
kind of fiow, namely the mean field EOS and the nucleon-nucleon scatterings, tends to be 
the same for a stiff, resp. soft EOS. In the soft case the scatterings dominate because of 
the larger compression reached but this is compensated by the less repulsive mean fields. 
In view of the failure to reproduce the observed size dependence it is possible to question 
this conclusion. The increased production of clusters in larger systems, correlated with 
increased radial fiow and stopping could well be a memory of the original compression 
(cooling by compression followed by expansion). In |3l] it was found that softer EOS 
favours cluster formation at freeze out. This interpretation is supported by the theoretical 
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Fig. 33. Radial flow energy in central collisions (6o < 0.25) as a function of system size. This ob- 
servable is deduced from the mass dependence (A=l,2,3) of average kinetic energies at (90 it 10)° 
cm.. The open (red) data points are from separate fits (see previous figures) varying both the 
offset and the slope per mass unit (i.e. the flow parameter). The full (black) data points are 
obtained requiring a common offset parameter (indicated in the panels as offset temperature Tq). 
The incident beam energies, varied over one order of magnitude, are indicated in the upper left 
corner of each panel. 
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Fig. 34. System size dependence of radial flow in central (6o < 0.25) collisions. The experimental 
data (black full triangles joined by straight lines) are compared with IQMD-SM predictions for 
OAA Gev (left) and 1.0^4 GeV (right) incident energies. 

work of ref . [61] . 



7.3 Heavier clusters 



For incident energies below OAAGeV we were able to obtain data for fragments up to 
nuclear charge 6-8, depending somewhat on the incident energy. Although such data, not 
separated by isotope, and more limited in acceptance (see Fig. |6]) have been published 
('PHASE r) and extensively discussed earlier by our Collaboration, |10|57f58f62ll63j . it is 
worthwhile to look again at these data in the light of the present more complete 'PHASE 
H' information and new insights. The relevant information is summarized in Fig. |35] where 
we show the nuclear charge dependence of kinetic energies in central collisions of Au+Au 
for five incident energies varying from 0.09A GeV to 0.4A GeV. An earlier evaluation of 
these data, showing similar features, has been presented in [63]. The results cannot be 
compared directly: in ref. [63j a coordinate system rotated into the flow axis was used and 
in-plane fragments were suppressed. The datasets presented in each panel of Fig. [351 least 
squares fitted by a linear function, hold for cm. polar angles of (5 — 45)° and (80 — 100)°, 
respectively. To allow a more direct approximate comparison with the mass-dependent 
data discussed before, we have multiplied the measured kinetic energies per nucleon by 
2Z and plotted them versus 2Z. Obviously the slopes of the fitted lines ('2Z method') are 
very different for the two polar angle ranges (in agreement with [15]), but they tend to 
gradually align as the energy is raised (see upper right panel in the figure), a confirmation 
of the global trend of increasing stopping in this range of incident energies (Fig. |2T|) . 
Increased stopping leads to a more isotropic distribution of kinetic energy, that was also 
visible for the light charged particles (Fig. |23l upper right panel showing the 0.4^4 GeV 
data). It is fair to remind at this stage, that while the data below 45° are well covered by 
our acceptance, the data around 90° are somewhat less reliable due to the necessity for 
a two-dimensional extrapolation (section H]). We are confident that the major conclusions 
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are not seriously affected as, below 0.2^4 GeV, our analysis was found to be consistent 
with the more complete INDRA data [33]. 



It is interesting to compare the apparent radial flow deduced by the '2Z method', Fig. 
with that deduced by the 'A3 method'. Fig. |3U This is shown in Fig. |36] for centrality 
bo < 0.15 (note that in order to match the centrality used in the analysis illustrated 
by Fig. [35l the centrality is higher than in Fig. [3T] which is constrained to 6o < 0.25 and 
therefore yields Ef values lower by 4 to 7 %). One sees the convergence of the two methods 
at 0.4^4 GeV where stopping is maximal, but at lower incident energies the A3 method 
yields higher values. We suggest to correlate this with the pronounced stopping hierarchy 
found at the lower energies (section[6]). Clearly then, one needs to correct the apparent mass 
dependences for this effect when converting them to radial flow. This requires transport 
codes that reproduce the observed clusterization and stopping hierarchy. 

At this stage we wish to make a comment to our flow analysis published in 1997, [40j. In 
the light of the present extensive information on incomplete stopping and hence strong 
deviations from isotropy, the earlier analyses assuming spheric expansion and neglecting 
partial transparency must be revised in the sense that the flow values found using the 
limited acceptance of the PHASE I of FOPI included a significant longitudinal (memory) 
component. Thus the analysis included some counter-streaming effects of the nucleons due 
to incomplete stopping (see Fig. |27]). As the method used in [IQ] was strictly constrained 
by energy conservation, the sum of radial and two-fluid counter flow turned out to be 
correct, but the interpretation of the mass or charge dependence of the kinetic energies at 
these forward angles as just radial (spheric) flow cannot be upheld. The data themselves 
remain valid, except for a subtle effect on the topology of the two-dimensional (mj vs 
rapidity i/z) spectra: applying the ERAT criterion on the apparatus limited in PHASE I 
to laboratory angles forward of 30° lead for the most central collisions (i.e. in the tails of 
the ERAT distributions) to an underestimation of the anisotropy as the selection criterion 
favoured events with a higher hit rate near the large angle apparatus limits, increasing the 
apparent isotropy. In the present work, with the significantly larger PHASE II acceptance, 
this effect, caused by event-by-event fluctuations, is minimized. 



7.4 Sytematics of equivalent temperatures 



In section [5] it was shown that the constrained rapidity distributions could be well repro- 
duced with a thermal ansatz, giving an equivalent temperature T^g. This represents an 
alternative way of systematizing mass dependences of phase space distributions. A sys- 
tematics of T^g is shown in the next four figures that close this section. We find a similarity 
with the systematics of average kinetic energies of fragments emitted around 90° : compare 
Fig. [3T]with Fig. [371 However, since the phase-space cuts are different and the stopping is 
incomplete (non-spherical topology), the deduced parameters are not substitutes for each 
other. While we have started with 90° kinetic energies primarily to join up to some of the 
earlier literature P5|42f45|l55] , the T^g characterizing the constrained rapidity distributions 
allow us to join up to higher energy heavy ion data, including some of our own work [15] . 
which are traditionally shown in (more appropriate) axially symmetric coordinate systems 
rather than spherical coordinate systems. The variation of Teg with particle type evidently 
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Fig. 35. Nuclear charge dependence of kinetic energies for various indicated beam energies in 
central collisions of Au+Au. The measured kinetic energies per nucleon are multiplied by twice 
the nuclear charge (i.e. by 2Z) and plotted versus IZ. Pink open squares: polar angles 5 to 45 
degrees, blue open triangles: polar angles 80 to 100 degrees. The straight lines are linear fits. The 
upper right panel shows the ratio of the deduced slopes (i.e. 'flow' energies) for the two types of 
polar angle intervals versus incident energy. 
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excludes a naive interpretation in terms of a common (equilibrium) temperature. How- 
ever, the observations illustrated in Figs. [H] to [16] go beyond just defining average kinetic 
energies since they allow a statement on the shapes of the distributions. 

In Fig. |37]we show data for 0.15 A GeV beam energy varying the system size. The same 
strong structures in the data for hydrogen and helium isotopes are visible and the A3 
method introduced earlier, joining the proton and the deuteron T^g with the average of 
the two mass three (^He and '^H) data points is seen to define an accurate slope that can 
be again taken to quantify radial flow. The lower right panel allows to assess the system 
size dependence. 

In Fig. [38] for data obtained at 1.0 A GeV, we confront the systematics for hydrogen 
isotopes (left panel) with that for pions [38] of both charges. Note the strikingly different 
system size dependence of nucleonic clusters and of pions: while the T^g of pions slowly (the 
data are plotted with an ordinate offset) decrease with system size, the contrary is true 
for the hydrogen isotopes. Again, as suggested earlier, this is correlated to the different 
consequences of compression (creation of pions) and decompression-cooling (reabsorption 
of pions, creation of clusters from nucleonic matter). 

Apparent temperatures for protons and deuterons in Ni+Ni collisions have already been 
published earlier by our Collaboration [IH]. If we interpolate the data in Fig. [38]to estimate 
the Ni+Ni values we obtain Tgg values exceeding those of |1H] by about 15%. It is unlikely 
that the difference is due to a different analysis method. One possible explanation is the 
different methods used to select the most central collisions. In [18] the multiplicity was 
maximized between the polar laboratory angles (7 — 30)°, while in the present work ERAT, 
determined using the full apparatus acceptance, was maximized. 



In Fig. [39] we show an observation typical for the SIS energy range: the presence of cold 
spectator matter in non-central collisions tends to cool locally the participant matter. 
We compare collisions of different systems at the same estimated value Apart of the 
participant matter. This effect can be studied in more detail by varying the angle relative 
to the reaction plane [28], but this will not be done here. 



45 





Ni+N 


O.ISAGeV 


bo<0.15 LCP 




60 


♦ He 




1 
♦ 


- 


50 


□ H 




^'^ 


- 


40 


- 


/-"""^ 


u 


- 




- ^^ 




♦ 


- 


30 


^ 






- 


20 


- 






^ 


10 


- 






- 



1 2 3 

mass (u) 





Au+Au0.15AGeV b(,<0.15 LCP 


80 


1 ' 1 ' 1 ' 
-♦He ♦ 
□ H /^ 




T ^^ 


60 




40 


_^J^ 


20 


- 








80 - 



60 



40 



20 



1 ' 


1 




♦ He 


♦ 




- QH 




^ 


- 


T,^ 


- 


-^ 


\ a 


♦ 


1 







12 3 4 

mass (u) 

0.15AGeVbn<0.15A3 



80 



60 



40 



20 



O Au+Au 
□ Xe+Csl 
ANi+Ni 




1 2 3 

mass (u) 



mass u 



Fig. 37. Equivalent transverse temperatures of hydrogen (open red squares) and helium isotopes 
(full blue diamonds) as function of their mass number in central Ni+Ni (upper left panel), 
Xe+Csl (upper right panel) and Au+Au collisions (lower left panel). A linear fit (A3 method) 
is also shown. The lower right panel allows to compare the three systems directly using the data 
points of the A3 method. 

Finally, we present in Fig. HO] for Au+Au at two different incident energies (0.4^4 and 
1.5 A GeV) a complete systematics of Teg for all five LCP varying the centrality. From the 
smooth trends the limits for 6o = can be deduced. Note the inversion of the hierarchy 
for the two He isotopes when passing from the higher to the lower energy. 



8 Chemistry 



The 47r reconstructions, explained in detail in section HI allow us to assess the integrated 
yields of all the ejectiles that constitute the bulk of the colliding system. For the most 
central collisions we present these yields in tabulated form in an appendix for the 25 
system-energies studied in this work. Except for the two lowest incident energies (O.OQA 
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Fig. 38. Equivalent transverse temperatures for symmetric central collisions at lA GeV beam 
energy. Left panel: data for hydrogen isotopes versus mass for three indicated systems. Right 
panel: data for pions versus system mass. Note the tendency (for the hydrogen isotopes) to 
converge to a common apparent offset 'temperature' T^ (left) and the systematic difference 
between the two charged pions, as well as their weak decreasing trend as a function of Agyg. 

and 0.12y4 GeV), we account for the total system charge with an accuracy of about 5% or 
better. Note that we include in this balance also the pions, [38] • Other created particles, 
such as kaons etc, are negligible on the 5% level (but are of course interesting for many other 
reasons). For the two lowest energies we miss significant contributions from heavy clusters 
beyond Z=8. The interested reader might consult INDRA-ALADIN works [33]l64|l65j for 
more complete distributions in the energy range at and below O.ISA GeV. 

Such data are often analysed using models assuming (local) chemical equilibrium with 
a unique 'chemical' temperature . An early, instructive presentation of the arguments in 
favour of a conjectured chemical equilibrium at freeze-out can be found in [66] . Even if one 
does not adopt the equilibrium assumption, (many of the observations presented in this 
work do not favour the equilibrium assumption) the degree of clusterization quantized by 
these yields can be used to constrain the outgoing (non-equilibrium) entropy [62|67] . This 
information, together with our extensive data on stopping, section [6], can be used to assess 
the viscosity to entropy ratio that is currently intensely investigated in connection with 
the quark-gluon phase thought to be copiously produced at the highest currently available 
energies [ 68|69|70j . The task of analysing our chemistry data using a modern thermal 
code that bridges the low-energy multifragmentation regime with the high energy particle- 
producing regime using a realistic approach for late decays is delayed to later publications. 
In [10] it was shown that a purely statistical interpretation of the non-collective part of 
the available energy using various statistical codes that included late evaporations was 
not successful as it strongly underestimated the yields of heavy clusters, such as oxygen. 
This conclusion is not touched by the now favoured interpretation of the collective energy 
as being to a sizeable part due to still two-fluid flow caused by incomplete stopping. In 
contrast to this failure one has to acknowledge a very reasonable rendering of the full 
nuclear charge distributions for Au on Au at 0.15 and 0.25^4 GeV when using the AMD 
code (Antisymmetrized Molecular Dynamics) [7T|72] . This transport code does not rely 
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Fig. 39. Equivalent transverse temperatures of protons (upper left panel), deuterons (upper 
right panel), tritons (lower right panel) and ^He as function of participant number Apart for 
collisions at E/A = 1.0 A GeV. The straight lines are linear fits. Open symbols and solid lines 
are for a variation of the impact parameter using Au+Au, full symbols and dashed lines are for a 
variation of the system size (Ca+Ca, Ru+Ru and Au+Au) at fixed reduced impact parameters 
(6o < 0.15). 

on the equilibrium assumption. 

Despite the different data analysis method, the yields listed in the appendix for Au+Au 
at 0.15, 0.25 and 0.4^4 GeV are fully compatible with those published earlier by our 
Collaboration |^0|f58] . The possibility to compare these yields with those published by other 
authors are limited: besides the need to match centralities, one is generally confronted with 
the lack of An estimates. As shown partially in |33], our data are nicely compatible with 
INDRA data in the low energy range where the measurements overlap and centralities 
have been carefully aligned. 

The situation on cluster production as it existed in the mid-eighties was described in 
Well suited for a direct comparison are yield ratios of light clusters (d, t, '^He, a) 
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Fig. 40. Centrality dependence of equivalent transverse temperatures for protons (circles), 
deuterons (squares), tritons (full triangles), ^He (diamonds) and ^He clusters in Au+Au col- 
lisions. Left panel: 0.4A GeV incident beam energy, right panel: 1.5A GeV. The smooth lines 
guide the eye and allow an estimate for bo = 0. 

Table 1 

Ratios of light cluster yields to proton yields in central collisions. Comparison of Plastic Ball 
data M\ with FOPI data. 





Au+Au 0.15A GeV 


Au+Au OAA GeV 


Nb(Ru)+Nb(Ru) 0.4>1 GeV 


ratio 


Plastic Bah 


FOPI 


Plastic Ball 


FOPI 


Plastic Bah 


FOPI 


d/p 


0.80 ± 0.08 


0.85 + 0.06 


0.65 + 0.07 


0.65 + 0.04 


0.50 + 0.05 


0.55 + 0.03 


t/p 


0.40 ± 0.04 


0.63 + 0.06 


0.40 + 0.04 


0.37 + 0.03 


0.28 + 0.04 


0.25 + 0.025 


3He/p 


0.50 ± 0.045 


0.35 + 0.04 


0.20 + 0.02 


0.18 + 0.02 


0.18 + 0.02 


0.17 + 0.017 


^He/p 


0.65 + 0.06 


0.69 + 0.07 


0.20 + 0.02 


0.22 + 0.02 


0.16 + 0.02 


0.16 + 0.016 



to protons measured [67] with the Plastic Ball spectrometer for the reactions Nb+Nb 
(which is close to one of our systems: Ru+Ru) and Au+Au for incident beam energies 
of 0.15, 0.25, 0.40 and 0.65 GeV/nucleon and were plotted by the authors versus the so- 
called participant proton multiplicity, Np, which was the total charge included in all the 
LCP cutting out some areas close to projectile and target rapidity to exclude 'spectators'. 
(Without these specifications the total charge would always be trivially constant in each 
collision due to charge conservation). Assuming Boltzmann like spectra, ratios were taken 
in limited regions of phase space where isotopes were identified and the momenta scaled 
with {1/mY'^. If we compare these ratios taken at the highest Np with corresponding ratios 
obtained from our Table for the most central collisions we get excellent agreement, for most 
cases. See the sample Table H] below for Au+Au at 0.15A, OAOA, and Nb+Nb/Ru+Ru at 
OAA GeV. 



The Plastic Ball data were read off the figures in [67] and the uncertainties somewhat arbi- 
trarily set at about 10% looking at data point straggling and assessing reading errors. The 
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agreement is not trivial as there are, besides the very different apparatus used, some differ- 
ences in the experimental procedure: In contrast to [67] we have obtained 47r reconstructed 
yields and no cuts on the data were done. Our procedure therefore does not involve the 
necessity to assume Boltzmann-like spectra, all with the same apparent temperature, nor, 
for these very central collisions, an arbitrary definition of 'spectator' contributions. Section 
|5] shows that such contributions would be difficult to identify in the rapidity distributions. 
The authors of ref. j6^, comparing with the predictions of a two parameter quantum sta- 
tistical code, QSM, [71], conjectured that the only serious discrepancy with the code (and 
as can be seen from Table H] also with our data), the t/p and ^He/p ratios for Au+Au at 
0.15A GeV, were due to deviations from the Boltzmann shapes. Our data fully support 
this conjecture: at the lower energies the spectra of ^He and ^H are very different, see also 
section [51 As mentioned before, we shall not repeat statistical model calculations in the 
present work although our data considerably extend the available information. Rather, it 
is desirable that such data be reproduced by the same transport model codes that are 
used to extract EOS information from the data: here the precise sorting with fixed cross 
sections of ERAT, an observable which is not directly connected with multiplicity, should 
be conceptually advantageous. 

Even without a detailed comparison with statistical or dynamical codes a number of 
interesting features can be deduced by just plotting some of the information enclosed in 
our chemistry data. We begin with some global features, shown in Fig. |4T1 As illustrated 
in the upper left panel, in the SIS energy range the percentage of single protons emitted 
in central Au+Au collisions starts at the lower end at roughly the 10% level to gradually 
increase to about 2/3 of the available charge at 1.5 A GeV. 

The fraction complementary to single protons, the clustered fraction, is plotted separately 
again in the lower left panel. Although it contains no additional information, since it must 
complement the proton fraction to 100%, the ordinate is now linear and the abscissa is 
now logarithmic, allowing, in this representation, to illustrate a linear trend above 0.2A 
GeV (see the straight line, that reproduces the data with an accuracy of 1.5%) that 
suggests by extrapolation that the clustered fraction stays above 10% all the way up to 
4:A GeV. This persistence of a significant probability to clusterize at freeze-out up to 
an available energy per nucleon more than two orders of magnitude higher than typical 
nucleonic binding energies is a remarkable signal of the local cooling process accompanying 
the fireball expansion and serves as a strong constraint on the associated entropy. 

In some of the literature, see for example [75j, particles with Z < 2 (LCP's) are supposed 
to be the ingredients of nuclear matter in the gas state. These are the particles most easily 
evaporated from a hot liquid and therefore forming the gaseous atmosphere above the 
liquid. If one accepts this simple view then the 'fifty-fifty' crossing liquid-gas takes place 
around 80A MeV beam energy. Since 'gas' and 'liquid' should add up to 100%, one can 
also read off the 'liquid' fraction from the figure, which we call heavy clusters (or 'droplets' 
|34] ) and which is called frequently IMF (intermediate mass fragments) in the literature. 
Above OAA GeV the liquid fraction is 5% or less. In the energy range considered here, fully 
central collisions do not lead any longer to fission-like (slightly less than half the projectile 
or target mass) or compound-nuclear-like (half or more than the total mass) remnants. 

The SIS energy regime is also the regime where particle creation, very dominantly pions. 
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Fig. 41. Some global characteristics of central Au+Au collisions as function of incident energy. 
Upper left panel: Reduced proton multiplicities and charge conversions by pion emission. Upper 
right panel: fraction (in %) of total energy emitted as pions. Lower left panel: % fraction of 
protons in clusters. Lower right panel: measured nuclear charge fraction (in %) emitted as LCP. 

varies from a 'perturbative' low level to a significant fraction expected to heavily influence 
the global dynamics. Thus, the energy fraction rises to the 20% level as seen in the upper 
right panel. The pertinent data were obtained from [38] . 



The non-perturbative nature of pions at the high end of this energy regime is also active 
for the isospin degree of freedom. According to the difference in tt" and vr^ production, 
[5S] . plotted also in the upper left panel, some eight protons are converted to neutrons at 
freeze out, modifying substantially the initial N jZ (from 1.494 to 1.264), a feature that in 
turn will modify the relative yields of the isotopes of hydrogen and helium. 

We will later compare yields of the isospin pair '^He and ^H. It is therefore of some interest 
to check if the yields are equal in an isospin symmetric system such as ^"Ca+^'^Ca. In 
this system only Coulomb effects could be invoked to predict different yields for the two 
isobaric clusters. Within experimental uncertainties the yields are equal in the investigated 
energy range (0.4 to 2 A GeV), see Fig. |l2l 
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We have already discussed in section |5] the non-perturbative character of clusterization. 
The high degree of clusterization seen in Fig. Hi] confirms this aspect. In Fig. US] we show 
excitation functions for the two He isotopes in central collisions of Au+Au. The two 
functions cross around 0.5A GeV. Clearly a simple perturbative coalescence model 
could not explain this behaviour. 



8.1 Mid-rapidity chemistry 



In the present energy range a clean separation between 'participating' (originally overlap- 
ping) matter and (peripheral) 'spectator' matter is not possible as, especially in the most 
central collisions, there are no readily identifiable 'spectators' moving undisturbed into an- 
gles close to the beam direction. Rather there is a superposition of various ejectiles having 
a more or less violent collision history as suggested by the observed stopping hierarchies 
(section [6]) . To minimize the effects of less violent and less equilibrated surface collisions, 
we concentrate our attention to (longitudinal) 'midrapidity chemistry', defined here by a 
cut \yzo\ < 0.5. As explained already in section [5l we remove some of the arbitrariness of 
the \yzo\ cut using the thermal model fits to the constrained rapidity distributions applying 
the same cut and then correcting for the cut-out part using the fitted model data. This 
gives us an estimate of mid-rapidity yields extended to full rapidity. 

In order to remove trivial size dependences we present ratios to proton yields, loosely 
writing e.g. d/p for the ratio of deuterons to protons. We shall see that even these 'size- 
reduced' observables (and hence the extracted effective thermal parameters) are definitely 
still system size dependent. Furthermore, using the isospin pairs ^H/^He and tt^/tt"*", we are 
able to separate size and isospin dependences: like the study of ground state masses, also 
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the variation of the reaction systems along the valley of stability to change the system size 
necessarily leads to strongly correlated simultaneous change of the isospin of the system. 
For 0.15y4, 0.4^4, l.OA and 1.5A GeV we show the results of this kind of analysis, see the 
next four figures, Figs. HUtoHZl and their captions. To each of these figures we add a panel 
showing the characteristics of the studied systems in terms of N/Z versus Asys- To better 
see the system size variations relative to the lightest system (Ni+Ni at the lowest energy, 
Ca+Ca at the higher energies), we also use double ratios. As a main conclusion from these 
rich data we can say that there is a significant influence of both system size and system 
composition on the outcoming yields of light clusters and of pions. When averaging over 
the isospin partners ^He/^H, resp. 7r~/7r+ we see (left lower panels in Figs. |45] to HTl) a 
regular evolution of the system-size dependence with the mass of the ejectile. In particular 
we confirm the weak dependence of total, or isospin averaged, pion multiplicities [38] . 

In contrast, in the upper two panels of the same figures the data for two isospin partners 
are seen to branch off as the N/Z is varied along with the size. In an isospin asymmetric 
{N > Z) medium like ^^^Au + ^^^Au (in contrast to ''^Ca + ""^Ca, Fig. |12]) more n-n 
than p-p collisions take place, hence more vr" than vr"'" are produced (see also our work 
ref. [38]). For similar combinatorics reasons more ^H than ^He are synthesized. The more 
interesting question is whether this is also influenced by density-dependent long-range 
isovector potentials. This question, so far, is not settled. For instance, for pion ratios, 
see the contradictig references [761I77II78] . Concerning ^He and ^H, to our knowledge no 
microscopic transport code has reproduced so far the special 'anomaly' features discussed 
earlier in connection with Figs. [T7] and [TS] in section 15.21 and Fig. ^U\ (right panel) in 
section [71 This makes also a convincing rendering of the experimental yield data difficult. 
Obviously more work is necessary. 

A final important point to make on the chemistry is the strong correlation of the yields, the 
transverse rapidity variances and the fragment specific stopping when varying the system 
sizes Zgyg. This effect is best illustrated when feeding from higher mass clusters is small so 
that one can envision the synthesis of new nuclei from a hot expanding, originally nucle- 
onic, soup. In view of the steep decrease of yields with fragment size, this condition could 
be safely assumed to be fulfilled at OAA GeV for Li [34j. As is shown in Fig. HHJthe same 
can be assumed to be true at l.OA GeV for ^He clusters and at 1.5 A GeV for mass three 
clusters. To remove superimposed isospin effects we have added the mass three H and He 
isotopes. Increased stopping (right panels) indicates increased compression. The increas- 
ing constrained transverse rapidity variances varxmO (middle panels) are interpreted as 
increasing radial flow developed thereafter in the expansion phase which is coupled to in- 
creased cooling ('droplet formation'). (We remind the reader that not only the variances, 
but also the shapes of these distributions were well described by a thermal shape assum- 
ing an equivalent temperature Teq, see section [5A] and flgs. 11-16.) The interpretation is 
supported by the remarkable, close-to-quadratic (the power is 1.7±0.1) system-size depen- 
dences of the mid-rapidity cluster multiplicities already stressed in ref. [51]. Kaon {K'^) 
production in the SIS energy range is also known [TH] to show stronger than linear size 
dependences associated to varying achieved compressions. The mechanisms are different 
however: in the kaon case there is increased production in the high density case which 
thereafter, in contrast to pion production Lll, is partially 'memorized' due to small K^ 



the insensitivity of pion yields to assumptions on the EOS and the system size is due, in this 
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Fig. 44. System size dependences for beam energies of 0.15A GeV. Lower left panel: yield ratios 
relative to protons for "^He (pink closed circles), ^He (red open diamonds), tritons (blue closed 
triangles) and deuterons (green open squares). Upper left panel: system size dependence of mass 
two, three and four yields relative to protons. The plotted double ratios are relative to the 
corresponding ratio in the Ni+Ni system {Agys = 116). The medium characteristics for the three 
studied systems in terms of N/Z and total mass Agyg are drawn in the lower right panel. 

reaction cross sections making the cooling phase 'inefficient'. For clusters we start 'from 
zero' at maximum achieved density, but increased compression leading to subsequent in- 
creased cooling creates more clusters. In this sense we suggest that this process, like K^ 
production, is sensitive to the EOS because it partially memorizes the initial compression 
which, in turn, is largest for softest EOS. 

Loosely speaking, we are thus witnessing the creation of nuclei in a mini- 'supernova'. 



interpretation, to a compensation of higher production at higher achieved maximum density by 
a stronger absorption-cooling when expanding 
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Fig. 45. System size dependences for beam energies of 0.4A GeV. Upper left panel: tjp and 
^He/p ratio, upper right panel: i:^ jp and 'n^ /p. The straight lines are linear fits excluding the 
Zr+Zr data points which are plotted with full symbols. The resulting (linear) 'isospin averaged' 
trends are also drawn. Lower left panel: double ratios (with respect to protons and relative to 
the Ca+Ca system, Agys = 80) for mass 4, 3, 2 fragments and pions (upper to lower curves). 
The Zr+Zr data points are not included here. The medium characteristics for the four studied 
systems in terms of N/Z and total mass Agys are drawn in the lower right panel. 



With increasing beam energy we catch earlier, hotter, stages forming smaller nuclei. For a 
theoretical treatment of the 0.4^4 GeV data in this connection see [61]. This work supports 
our interpretation and shows that understanding this 'fast' synthesis of nuclei is not outside 
our present theoretical capabilities, although there is room for improvement. 
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Fig. 46. System size dependences for beam energies of l.OA GeV. Upper left panel: t/p and 
^He/p ratio, upper right panel: vr'^/p and ir^ /p. The straight lines are linear fits guiding the eye. 
The resulting (linear) 'isospin averaged' trends are also drawn. Lower left panel: double ratios 
(with respect to protons and relative to the Ca+Ca system, Agyg = 80) for mass 4, 3, 2 fragments 
and pions (upper to lower curves). The medium characteristics for the three studied systems in 
terms of N/Z and total mass Agyg are drawn in the lower right panel. 

9 Summary 



We summarize here the results presented in sections |5]to[8l In a study involving 25 system- 
energies we have been able to establish a vast systematics of very central heavy ion col- 
lisions in the SIS energy range. We started characterizing the data for identified ejectiles 
by rapidity distributions both in the longitudinal and the transverse directions relative 
to the beam axis. As a very general rule we found that the variance of the longitudinal 
distribution was always broader than the variance of the transverse direction. This effect 
was most pronounced at the highest incident energy (Figs. [T H[T2|) . 
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Fig. 47. System size dependences for beam energies of 1.5A GeV. Upper left panel: tjp and 
^He/p ratio, upper right panel: i:^ jp and 'n^ /p. The straight lines are linear fits guiding the eye. 
The resulting (linear) 'isospin averaged' trends are also drawn. Lower left panel: double ratios 
(with respect to protons and relative to the Ca+Ca system, Asys = 80) for mass 4, 3, 2 fragments 
and pions (upper to lower curves). Data for the Zr+Zr system are drawn as full symbols and 
are not included in the fits. The medium characteristics for the four studied systems in terms of 
N/Z and total mass Agys are drawn in the lower right panel. 



We also introduced a transverse distribution constrained by a cut \yzo\ < 0.5 on the scaled 
longitudinal rapidity (the definitions were given in section [3]) and found that its shape, 
plotted on a linear scale, could always be very well described by a thermal distribution 
defining thus an equivalent 'temperature' T^q that was characteristic for each particle, 
incident energy and centrality. This equivalent temperature was found in general to rise 
with the particle's mass. Fig. [131 an effect that could be qualitatively reproduced by 
IQMD simulations. Fig. [TH although the shapes of the simulated distributions looked less 
perfectly thermal and the Teq followed the experimental data only approximately. 
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Fig. 48. System size dependences of symmetric central collisions at incident beam energies of 
l.OA GeV (lower panels, ^He fragments) and 1.5A GeV (upper panels, average of ^He and ^H 
fragments). All straight lines are linear fits. Left panels: Multiplicities (reduced to 100 protons). 
The data marked 47r with open square symbols are 47r integrated, while the closed symbols 
represent mid-rapidity yields (constrained to ly^ol < 0.5). Middle panels: Variance varxmO of the 
scaled constrained transverse rapidity distributions (dn/di/xmo)- Right panels: stopping varxz. 

At 1.5 A GeV beam energy the rise of T^q with mass could also be followed with the 
He isotopes and as a consequence of this behaviour it could be shown that the spectra 
of heavier particles could not simply be derived by power law from the proton spectra, 
Fig. [15] somewhat in contradiction to messages published in ref . [43j . 

The Teq were shown to rise significantly with centrality. Fig. [161 At lower incident energies 
a remarkable deviation of the mass trend of Tgq was observed: the '^He value tended to be 



lower than the ^He value. This 'He anomaly' [S] and large differences between ^H and '^He 
spectra were tentatively associated with clusterization phenomena, Figs.[T71[IHl illustrating 
the fact that a proper rendering of cluster probabilities was essential if simulations were 
to be used for extracting basic information from the data. 



The systematic comparison of transverse and longitudinal rapidity distributions. Fig. [191 
turned out to be a rich source of information, as one observes remarkable variations with 
incident energy and system size. Although the shapes of the distributions showed features 
that cannot be fully described by one parameter, we proposed to summarize the differ- 
ential information by a stopping observable varxz, the ratio of variances of transverse to 
longitudinal rapidity distributions. Stopping was always found to have values below one, 
a fact that we interpreted as partial transparency or incomplete stopping on account of 
the system size dependences. The hydrodynamic rebound expected in the full stopping 
scenario could also be definitely excluded by use of isospin tracer methods introduced ear- 
lier [52] by our Collaboration: an increase of transparency at the high end of the studied 



energy range, Fig. [201 could be unambiguously confirmed. 

We were able to confront the global stopping observable [32] with simulations varying 
the stiffness of the nuclear EOS. The most striking feature of the excitation function for 
the Au+Au system, namely the existence of a pronounced maximum for E/A = 0.2A to 
0.8A GeV could be qualitatively reproduced, but the strong descent at higher energies 
to decreased stopping was underestimated. Fig. |2T1 The predictions for varxz were found 
to be sufficiently sensitive to the EOS to exceed the experimental uncertainty: this adds 
another constraint from central collisions to efforts to fix the EOS from heavy ion data. 

Moving beyond global observables to individual stopping data associated to specific iden- 
tified particles, we found stopping to be essentially a collective phenomenon. Fig. |221 but 
a closer inspection also established that there was a hierarchy of stopping suggesting that 
heavier clusters were in part connected on the average with a less violent collision his- 
tory. This hierarchy was weakest in the maximum of the stopping excitation function. 
Fig. [231 but otherwise followed as a function of isotope mass, incident energy and system 
size a rather complex behaviour reflecting the rich data presented earlier in Fig. [19] and 
representing a challenge to future simulations. Fig. [2l] Surprisingly, we could not find a 
significant dependence of stopping on isospin at 0.4^4 and 1.5A GeV, Fig. [25] while the 
expected strong dependence on centrality. Fig. [26] is a warning that centrality must be 
well matched to experiment when simulations are done. 

For the lower energies we have confirmed that the stopping hierarchy is omnipresent also 
for all heavier clusters up to Z = 8. In an attempt to join up to energies higher than were 
studied here we showed the strong trend for increasing transparency of leading protons as 
the energy is raised to AGS and SPS energies. 

Passing from the compression stage, characterized by the stopping observable, we studied 
the expansion stage, characterized by radial fiow, which is generally inferred by mass 
dependences of kinetic energies or equivalent temperatures. We found that average kinetic 
energies of identified particles emitted at 90° follow extremely regular patterns as a function 
of incident energy. Fig. [29] the details of which, like the ^H-^He difference, can only be 
understood if the simulation predicts clusterization correctly. Figs. [30l[3l] If we combine 
the data for protons, deuterons and the average of the two mass three isotopes ('A3 
method'), we can establish a radial flow systematics for Au+Au extending from 0.09A to 
1.5 A GeV which is in fair agreement in the overlapping energy range with [15]. We flnd 
a strong system size dependence of radial flow. Fig. [331 which is however not reproduced 
quantitatively by IQMD simulations. Fig. [HH In agreement with [45] we flnd that the mass 
dependences of kinetic energies are strongly dependent on the emission angles. Fig. [35] 
leading to a necessary reassessment of the origin of the 'flow' deduced in our earlier work on 
heavy clusters [IQ] . In general the observed partial transparencies and stopping hierarchies 
call for a more qualifled interpretation of the data, as illustrated in Fig. 



Using the equivalent temperatures, based on constrained transverse rapidity distributions, 
we observe the same structural characteristics if we plot Tgg versus mass for the flve LCP, 
Fig. [371 The system size evolution of T^q is quite different for nucleonic clusters than for 
pions. Fig. [38] reflecting the different mechanisms of creation/annihilation. The presence 
of 'spectator' matter in non-central collisions influences the resulting Teq, Fig. 
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A systematics of T^g for all five LCP versus impact parameter is presented in Fig. HO] and 
allows to extrapolate to 60 = 0. 

The final state of the system after expansion-cooling is described by its chemical compo- 
sition which we have tabulated in the appendix for the twenty-five system energies. For 
energies below 0.2A GeV we are consistent with INDRA-ALADIN data [33|64|I65] . while 
quite satisfactory agreement with the pioneering PLASTIC BALL experiments [67] could 
be demonstrated, see Tabled for cases of overlap. Throughout the SIS-BEVALAC energy 
range clusterization is an important feature, Fig. HH that cannot be treated by perturba- 
tive coalescence models and influences significantly the spectra and flow of all particles: 
hence conclusions from the data with the help of transport codes must handle correctly 
this aspect of the reactions. At the same time, particle creation (pions) in this energy 
range rises gradually from a perturbative level to as much as 20% in terms of the energy 
fraction at 1.5 A GeV (many more details, connected with pions and obtained with the 
FOPI setup, were published in [38] )• At the lower energy end a particles play a special 
role: below 0.5A GeV more a particles are emitted than ^He ejectiles. Fig. |131 and a par- 
ticle condensation is probably needed to explain in detail the very different ^He and ^H 
spectra as suggested by our discussion of FigsHT] and [181 There is a notable system-size 
and isospin dependence of the chemical 'temperatures' that we evidenced using isospin 
pairs like '^He-'^H and 7i~^-7i~ in various systems. Figs. H^ll471 



Finally, our yield data suggest that there is at least a partial memory of the initial degree of 
compression, and hence the stiffness of the EOS: more compression leads to more efficient 
subsequent cooling enhancing the degree of clusterization. This is suggested by the close 
to quadratic dependence on system size and the conspicuous correlations between cluster 
yields, radial flow and stopping evidenced in the last figure of this paper and in 
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Appendix: Yield tables 

Multiplicities and charge balance for Au+Au at E/A = 0.09 GeV and Bq < 0.15. 

Z=l 40.1 ±2.1 p 17.0 ±0.9 Z=2 23.2 ± 1.7 ^He 
d 13.3 ±1.0 ^He 

t 9.7±1.0 
Li 4.1 ±0.4 Be 1.77 ±0.18 
B 1.79 ±0.18 C 1.43 ±0.15 



Charge balance LCP: 86.5 or 54.8% 

Charge in HC Z3-Z6 (measured): 36.8 ± 3.7 or (23.3 ± 2.4)% 

Charge measured: 123.3 or 78.0% 

Charge in HC: 158-86.5=71.5 or 45.3% 



Multiplicities and charge balance for Au±Au at E/A = 0.12 GeV and Bq < 0.15. 

Z=l 51.5 ±2.6 p 21.5 ±1.1 Z=2 24.2 ±1.7 ^He 8.0 ± 0.8 

d 17.6 ±1.2 ^He 16.2 ± 1.6 

t 12.2 ±1.2 

Li 4.7 ±0.5 Be 1.75 ±0.18 

B 1.62 ±0.16 C 1.16 ±0.12 

N 0.71 ±0.07 O 0.37 ±0.04 



Charge balance LCP: 99.7 or 63.1% 

Charge in HC Z3-Z8 (measured): 43.9 ± 4.4 or (27.8 ± 2.8)% 

Charge measured: 143.6 or 90.9% 

Charge in HC Z9-Z14 (estimated): 8.2 ± 3.0 or (5.2 ± 1.9)% 

Total balance: 96.1% 

Adopted IMF charge: 55.2 ± 5 or 34.9% 



Multiplicities and charge balance for Ni±Ni at E/A = 0.15 GeV and Bq < 0.15. 

Z=l 27.3 ±1.4 p 15.1 ±0.8 Z=2 10.8 ± 0.8 ^He 4.06 ± 0.41 

d 8.4 ±0.6 ^He 6.71 ± 0.68 

t 3.74 ±0.4 

Li 1.24 ±0.13 Be 0.42 ± 0.04 

B 0.25 ±0.03 C 0.18 ±0.02 
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Charge balance LCP: 48.84 or 87.2% 
Charge in HC: 7.69 or (13.7 ± 1.4)% 
Total balance: 56.53 or 100.9% 



Multiplicities and charge balance for Xe+Csl at E/A = 0.15 GeV and b^ < 0.15. 

Z=l 46.0 ±2.3 p 20.1 ±1.0 Z=2 18.8 ± 1.3 ^He 5.9 ± 0.6 

d 15.6 ±1.1 ^He 12.9 ±1.3 

t 10.3 ±1.0 

Li 2.77 ±0.28 Be 0.89 ± 0.09 

B 0.67 ±0.07 C 0.41 ±0.04 

N 0.19 ±0.02 O 0.10 ±0.01 



Charge balance LCP: 83.53 or 77.3% 
Charge in HC: 19.83 or (18.4 ± 1.8)% 
Total balance: 103.36 or 95.7% 



Multiplicities and charge balance for Au±Au at E/A = 0.15 GeV and Bq < 0.15. 

Z=l 59.9 ±3.0 p 24.1 ±1.2 Z=2 25.0 ±1.8 ^He 8.4 ± 0.9 

d 20.4 ±1.4 ^He 16.6 ± 1.7 

t 15.1 ±1.5 

Li 5.0 ±0.5 Be 1.69 ±0.17 

B 1.44 ±0.15 C 0.90 ±0.09 

N 0.50 ±0.05 O 0.26 ±0.03 



Charge balance LCP: 109.9 or 69.6% 

Charge in HC Z3-Z8 (measured): 40.0 ± 4 or (25.3 ± 2.5)% 

Charge measured: 149.9 or 94.9% 

Charge in HC Z9-Z14 (estimated): 4.4 ± 1.5 or (2.8 ± 1)% 

Total balance: 97.7% 



Multiplicities and charge balance for Ni±Ni at E/A = 0.25 GeV and bo < 0.15. 

Z=l 34.9 ±1.8 p 19.2 ±1.0 Z=2 9.0 ± 0.7 ''^He 3.24 ± 0.33 
d 10.5 ±0.8 ^He 5.79 ± 0.58 

t 5.1 ±0.5 
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Li 0.91 ±0.09 Be 0.26 ± 0.03 
B 0.10 ±0.01 

Charge balance LCP: 52.93 or 94.5% 
Charge in HC: 4.28 or (7.6 ± 0.8)% 
Total balance: 57.21 or 102.2% 

Multiplicities and charge balance for Xe±CsI at E/A = 0.25 GeV and bo < 0.15. 

Z=l 59.6 ±3.0 p 26.3 ±1.4 Z=2 16.7 ±1.2 ^He 6.0 ± 0.6 

d 19.3 ±1.4 ^He 10.8 ± 1.1 

t 14.1 ±1.5 

Li 2.33 ±0.23 Be 0.60 ± 0.06 

B 0.29 ±0.03 C 0.14 ±0.02 

N 0.05 ±0.01 



Charge balance LCP: 92.99 or 86.1% 
Charge in HC: 12.65 or (11.7 ± 1.2)% 
Total balance: 105.64 or 97.8% 



Multiplicities and charge balance for Au±Au at E/A = 0.25 GeV and 6o < 0.15. 

Z=l 79.9 ±4.0 p 33.7 ±1.7 Z=2 24.8 ±1.7 ^He 8.8 ± 0.9 

d 26.7 ±1.9 ^He 16.0 ± 1.6 

t 19.5 ±2.0 

Li 4.5 ±0.4 Be 1.32 ±0.13 

B 0.67 ±0.07 C 0.35 ±0.04 

N 0.13 ±0.02 



Charge balance LCP: 82.1% 

Charge in HC: 26.0 or (16.45 ± 1.7)% 

measured HC charge (Z3-Z7): 25.08, estim (Z8-Z10): 0.91 

Total balance: 98.5% 



Multiplicities and charge balance for Au±Au at E/A = 0.40 GeV and feg < 0.15. 
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Z=l 106.0 ±4.0 p 52.9 ±2.7 Z=2 21.3 ± 1.9 ^He 9.4 ± 0.9 
d 34.2 ±2.1 '^B.e 11.9 ±1.1 

t 19.4 ±1.8 
7r+ 0.95 ±0.08 Tr" 2.80 ±0.14 

Li 3.5 ±0.4 Be 0.84 ± 0.09 
B 0.27 ±0.03 C 0.096 ±0.01 



Charge balance LCP: 148.6 or 94.1% 
Charge balance LCP±pions: 146.8 or 92.9% 
Charge in HC (Z3-Z6): 15.77 ± 1.6 or (10.0 ± 1.0)% 
Total balance: 102.9% 



Multiplicities and charge balance for Au±Au at E/A = 0.6 GeV and Bq < 0.15. 
Z=l 121.9 ±6.1 p 68.5 ±3.5 Z=2 16.4 ± 1.5 ^He 9.1 ± 0.9 
d 36.6 ±2.2 ^He 7.4 ± 0.7 

t 16.8 ±1.5 
7r+ 3.14 ±0.25 tt" 7.08 ± 0.36 



Charge balance LCP±pions: 95.4% 
Estimate IMF: (6 ± 2)% 
Total balance: 101.4% 



Multiplicities and charge balance for Au±Au at E/A = 0.8 GeV and bo < 0.15. 
Z=l 132.51 ±6.6 p 79.4 ±4.0 Z=2 13.5 ± 1.2 ^He 8.2 ± 0.8 

d 37.7 ±2.3 ^He 5.2 ± 0.5 

t 15.4 ±1.6 
7r+ 6.26 ±0.50 tt" 12.49 ± 0.63 



Charge balance LCP±pions: 153.3 or 97.0% 
Estimate IMF: (4 ±2)% 
Total balance: 101.0% 



Multiplicities and charge balance for Au±Au at E/A = 1.0 GeV and bo < 0.15. 
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Z=l 140.9 ±7.0 p 91.5 ±4.0 Z=2 11.0 ±1.2 ^He 7.3 ± 0.7 

d 38.8 ±2.3 ^He 3.7 ±0.4 

t 13.9 ±1.3 

7r+ 9.87 ±0.79 tt" 18.26 ± 0.91 

Charge balance LCP±pions: 154.5 or 97.8% 
IMF order 2% (< 4%) 
Total balance: 99.8% 

Multiplicities and charge balance for Au±Au at E/A = 1.2 GeV and Bq < 0.15. 



Z=l 


148.2 ± 7.0 


P 


99.3 ±4.0 


Z=2 8.74 ±1.3 


3He 


6.1 ±0.9 






d 


36.8 ±2.3 




^He 


2.6 ±0.4 






t 


12.1 ±1.3 








7r+ 


14.13 ±1.13 


TT^ 


24.37 ±1.22 









Charge balance LCP±pions: 155.5 or 98.4% 

IMF neglected 

Total balance: 98.4% (IMF neglected) 



Multiplicities and charge balance for Au±Au at E/A = 1.5 GeV and bo < 0.15. 

111.3 ±5.6 Z=2 6.76 ±1.3 ^He 5.1 ± 0.8 
36.0 ±2.2 ^He 1.7 ±0.3 

11.0±1.1 
33.34 ±1.67 



L=l 


158.4 ±7.0 


P 
d 

t 


7r+ 


21.16 ±1.70 


TT" 



Charge balance LCP±pions: 159.7 or 101.1% 

IMF neglected 

Total balance: 101.3% (IMF neglected) 



Multiplicities and charge balance for Ru±Ru at E/A = 0.4 GeV and bo < 0.15. 
Z=l 59.8 ±3.0 p 33.2 ±1.7 Z=2 11.0 ± 1.0 ^He 5.6 ± 0.5 

18.3 ±1.1 ^He 5.4 ±0.5 

8.3 ±0.8 
7r+ 0.70 ±0.06 n- 1.07 ±0.06 

Li 1.28 ±0.13 Be 0.28 ± 0.03 
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59.8 ±3.0 


P 




d 




t 


0.70 ±0.06 


71' 



Charge balance LCP+pions: 81.37 or 92.5% 

Total charge IMF : 5.5 or 6.2% (estimate 0.5 for Z > A) 

Total balance: 92.5 + 6.2 = 98.7% 



Multiplicities and charge balance for Zr+Zr at E/A = 0.4 GeV and Bq < 0.15. 
Z=l 56.0 ±2.8 p 28.8 ±1.5 Z=2 10.2 ± 1.0 ^He 4.8 ± 0.4 

18.3 ±1.1 ^He 5.4 ±0.5 

8.9 ±0.8 
7r+ 0.56 ±0.05 vr" 1.42 ±0.07 

Li 1.22 ±0.13 Be 0.29 ± 0.03 



56.0 ±2.8 


P 




d 




t 


0.56 ±0.05 


71' 



Charge balance LCP±pions: 75.49 or 94.4% 

Total charge IMF : 5.3 or 6.6% (estimate 0.5 for Z > A) 

Total balance: 94.4 ± 6.6 = 101.0% 



Multiplicities and charge balance for Ru±Ru at E/A = 1.0 GeV and bo < 0.15. 
Z=l 79.1 ±4.0 p 53.9 ±2.7 Z=2 5.17 ±0.47 ^He 3.68 ± 0.33 
d 19.7 ±1.2 ^He 1.49 ±0.14 

t 5.5 ±0.5 
71+ 6.05 ±0.49 TV- 8.09 ±0.41 



Charge balance LCP±pions: 87.38 or 99.3% 
IMF neglected 



Multiplicities and charge balance for Ru±Ru at E/A = 1.5 GeV and bo < 0.15. 
Z=l 87.03 ±4.4 p 64.5 ±3.2 Z=2 3.32 ± 0.3 ^He 2.65 ± 0.24 

d 18.4 ±1.1 ^He 0.67 ±0.06 

t 4.1 ±0.4 
7r+ 12.49 ±1.00 n- 15.63 ±0.78 

Charge balance LCP±pions: 90.52 or 102.9% 
IMF neglected 



Multiplicities and charge balance for Zr±Zr at E/A =1.5 GeV and bo < 0.15. 
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Z=l 


30.5 ±1.6 


P 
d 

t 


7r+ 


0.45 ± 0.04 


71' 



Z=l 84.3 ±4.2 p 61.3 ±3.1 Z=2 3.26 ± 0.32 ^He 2.61 ± 0.24 

d 18.6 ±1.1 ^He 0.65 ±0.06 

t 4.3 ±0.4 

7r+ 11.45 ±0.92 tt" 17.19 ±0.86 

Charge balance LCP±pions: 85.09 or 106.4% 
IMF neglected 

Multiplicities and charge balance for Ca±Ca at E/A = 0.4 GeV and Bq < 0.15. 

19.1 ±1.0 Z=2 4.9 ±0.5 ^He 2.8 ± 0.3 

8.6 ±0.5 ^He 2.1 ±0.2 

2.7 ±0.3 
0.47 ±0.03 

Charge balance LCP±pions: 40.29 or 100.7% 

Z3 0.34 ±0.04 

Charge balance: 40.29±1. 02=41.31 or 103.3% 

Multiplicities and charge balance for Ca±Ca at E/A = 0.6 GeV and Bq < 0.15. 
Z=l 33.2 ±1.7 p 22.9 ±1.2 Z=2 3.7 ±0.4 ^He 2.6 ± 0.3 
d 8.1 ±0.5 ^He 1.1 ±0.1 

t 2.2 ±0.2 
7r+ 1.13 ±0.09 n- 1.28 ± 0.07 

Charge balance LCP±pions: 40.54 or 101.3% 
IMF neglected 

Multiplicities and charge balance for Ca±Ca at E/A = 0.8 GeV and Bq < 0.15. 

25.3 ±1.3 Z=2 2.6 ±0.2 ^He 1.85 ±0.23 

^He 0.75 ±0.09 



Z=l 


35.2 ±1.8 


P 


25.3 ±1.3 






d 


8.0 ±0.5 






t 


1.9 ±0.2 


7r+ 


2.11 ±0.17 


TT^ 


2.28 ±0.11 



Charge balance LCP+pions: 40.24 or 100.6% 

Multiplicities and charge balance for Ca±Ca at E/A = 1.0 GeV and 69 < 0.15. 
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Z=l 36.5 ±1.8 p 26.5 ±1.4 Z=2 2.1 ±0.2 ^He 1.65 ±0.15 

d 7.8 ±0.5 ^He 0.45 ± 0.04 

t 1.52 ±0.14 

7r+ 3.16 ±0.26 tt" 3.33 ±0.17 

Charge balance LCP±pions: 40.57 or 101.4% 

Multiplicities and charge balance for Ca+Ca at E/A = 1.5 GeV and Bq < 0.15. 
Z=l 38.6 ±1.8 p 30.4 ±1.6 Z=2 1.10 ±0.16 ^He 0.93 ±0.16 
d 7.2 ±0.5 ^He 0.17 ±0.03 

t 1.22 ±0.15 
71+ 5.62 ±0.45 TT" 5.78 ±0.30 

Charge balance LCP±pions: 40.64 or 101.6% 

Multiplicities and charge balance for Ca±Ca at E/A = 1.93 GeV and Bq < 0.15. 
Z=l 39.6 ±1.8 p 33.7 ±1.7 Z=2 0.82 ± 0.20 ^He 0.70 ±0.14 
d 6.8 ±0.5 ^He 0.13 ±0.03 

t 1.11±0.15 
7r+ 7.82 ±0.63 tt" 7.86 ± 0.40 

Charge balance LCP±pions: 41.22 or 103.1% 
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